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ABSTRACT 


This work examines hyperthermia as a cancer treatment 
and describes the development of a new hyperthermia 
applicator. The r.f. applicator produces more heat deep 
Within “a =Svbiect than on or near the surface sccording to 


the results obtained. 


Experimental temperature distribution results in an 
agar cylinder were found to be in close agreement with 
computer-generated data. The effects of altering various 
applicator parameters were explored uSing the computer 
models. From this work an optimized applicator was found 
which should produce more heat deep within a dielectric 


than on the surface. 


The hyperthermia applicator consists of parallel plate 
electrodes the same width as the heated medium (agar). 
Radio: Epequency. energy ati )3./75 MHz, or dower, is used 
because OT its substantial penetration depth. A 
temperature controlled saline bath is utilized to cool the 
surface of the irradiated dielectric and improve energy 
coupling. Finally the electrodes are made to orbit the 
irradiated medium to reduce superficial energy exposure 


without affecting deep exposure within the medium. 
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Live animal experiments were performed but proved to 


be inconclusive. 
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Calculated temperature distribution using the 
Same parameterS as in Figure 4.5 except that the 
electrodes revolve completely around the phantom. 


eeoeee#eege oeeee#ee oeee3sc5eeeee c PAVE ATRL Ra here Cee e tere acca to oeesieibee 


Calculated temperature distribution using the 
Same parameters as in Figure 4.9 (asymmetric 
electrodes) except that the electrodes revolve 


completely around the phantom... ..cesceeseeeee 146 


Xvi 


a : 
at yresepud hedsivsie) 
7 = © dgasqees . 2 . . on: = phi thi 


: Ry 


_ 


: RH 0 coy st o whens So 4456 fovcigeanté ga:0T 


— 
uae gates eoisudiaealhs siese5%5>>° Satotunt as 


: 7 
we} iene “6E,h. SHQEy wi 28 ireitranaq sake 
i) 


De ee ed ae Pr. s°Sa sate l33209 
1! 
| = gnive setzsudeiee saud: sipet bese, use. 
' ‘as tats aqwive $6.5 eure fi 6 cuwresases iptiz 


.) ee ie wate a eee Ww ue . = het Ae 306.2 


ats gia seas sis =eg@s: = 6&ee at eS 
wd aia cal? cheese D2, @idgiT om: 8) 2: *Seraag; see 
’ gebtaeta) OAS Onin BAe gees se 4°71 Seng dss 


Peery a. ° e ee 


7 | 

mg gate sod seh y sai pabtovedast, ibs istaolas 
. o) bt amped i. 26 .S1e2¢8e7Ad) ends 
hdily selwe (>ober2sefs 


ape ad + onion! watz Gniic-s «lade iggaa>)*' 


shin sd getae 


Ae 


Calculated temperature distribution using the 
same parameters aS in Figure 4.45 (9.7 cm wide 
Stationary symmetric electrodes) except that the 


electrodes rotate completely around the phantom.. 
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INTRODUCTION 


A HISTORY OF HYPERTHERMIA 


Heat has long been used by man in the treatment of 
disease. Heat in the form of naturally occurring fever or 
inflammation was naturally used to help control many 
infections before man learned to produce it. Hippocrates 
once said,"Give me the power to produce fever and I will 
cure all disease" [cited in 54]. Many afflictions have 
been, or are, treated with heat. The list includes acne 
loole arthritis (66), cancer a7], carbuncles and furuncles 
(iol eeecivabetes  [oi1ted sineco)], hemorrhoids: (25) scyonhy.tis 


[11], and possibly the common cold [43]. 


Asse long, Pago -as = s0U0T B.C... heat im the storm ofa 
smoldering stick was used to treat tumors [10]. In the 
1890's Coley showed that pyrogens (Coley's Toxins) could 
produce regression or even eNre neoplastic lesions [17]. 
Unfortunately he could not tell whether heat or some other 
immune response was responsible for the tumor regressions 
P70 ne Around the same time D'Arsonval and Tesla suggested 
and examined the use of high frequency (h.f.) fields (1-10 
MHz) in medicine [66]. Such effects as analgesia, vascular 
dilatation, and lower blood pressure were reported by 
D'Arsonval in 1893 [6]. “A colleague of D’Arsonval’s, Paul 
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purpose of treating tumors [cited in 66]. In 1906 de 
Keating-Hart used h.f.fields in the form of an electrical 


discharge to destroy tumors [40]. 


Interest in hyperthermia aS a cancer treatment was 
increased in 1918 by Rohdenburg. He published a literature 
review showing that half the recorded spontaneous tumor 
regressions were preceded by high fever, infection, or heat 
application ii]. pThiceindicated that some icancers could) be 
treated by any method which would increase the tumor 


temperature sufficiently. 


in 1939 Denier explored the use of microwave (3.75 
MHZ) @idiathermy 2m conjunction with » x-rays, ©as a tumor 
therapy [clted in 40]. yim saddition, the ~possibility® of 
selective heating as a result of variations in frequency 
and dielectric properties of tumors and organs had _ been 
eonsidqerea “by 1941 [cited in 966). Up to 1955 several 
papers were published dealing with the effects of heat and 
heat plus ionizing radiation on cancers. Interest in 
hyperthermia then waned until the early 1960's, when Crile 
published results showing that heat and especially heat 
plus ionizing radiation were effective in treating tumors 
(alice More recently, work by Muckle and Dickson [20, 46], 
Pettigrew [51], LeVeen [38, 39], Cavaliere [12], and many 
others, have led to the present interest in hyperthermia as 


a cancer therapy. 
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To date numerous papers have been published on the use 
and effect of heat on cancer tumors. The overwhelming 
conclusion is that hyperthermia, particularly when used in 
conjunction with chemotherapy or ionizing radiation, has a 
detrimental effect on tumors while sparing healthy tissue 
fase S524¢ There are also accounts in the literature 
claiming hyperthermia can act aS an analgesic and cause 
regression or eradication of tumors which were unresponsive 
to conventional therapy [50, 37]. The evidence shows’ that 
hyperthermia 1S a promising tool for cancer therapy and 


must by=anvestigated further. 
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THE EFFECT OF HYPERTHERMIA ON TUMORS 


At “a biochemical level, Ludgate. et aj [41] have 
described the effect of hyperthermia on tumor cells. Heat 
destroys sensitive cells by scrambling molecular bonds. 
This leads to changes in the cytoplasm and DNA damage. The 
DNA damage is not repaired because of heat induced changes 
in repair enzymes and their production. The cell walls 
become leaky, releasing ions. Aerobic glycolysis gives way 
tOmranaerobre: (toliycolysis Tesultang Win aweivower “pH “and 


Subsequent lysosomal activity. 


Tumors often consist of a well-vascularized periphery 
Surrounding a poorly vascularized, semi-necrotic, centre 
[38,54]. | Thescentral portion of the ‘tumor tends to be most 
vulnerable to heat therapy and resistant to ionizing 
radiatron[47,.48).8 Thercpposite as true for the: well 
Vascularized § pontion Wolm@ithe? tumor...) Growing cell's an the 
S-phase of mitosis have been found most vulnerable to heat 
WhilemnesiStant: toMmonigingemadratwonm {1.4 7)ee/2 lea) Hence: 
a fortunate synergism between hyperthermia and ionizing 
hadtation plexists 611468 5% Some researchers have found 
that certain chemotherapeutic agents appear to increase the 
Sensitivity of cells) to heat and vice versa {45)067]. In 
addition, there may be a tumoricidal immunological response 


associated with hyperthermia [65, 38, see p. 15]. 
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Investigations by LeVeen et a][38] have shown that 
growing tumors will displace normal vasculature and produce 
their own blood vessels through angiogenesis. This process 
occurs primarily at the periphery of the tumor [54]. These 
new vessels tend to be a mass of small capillaries which 
offer a high resistance to heat dissipating and nutrient 
beering “blood flow. "As a result, the centre of the tumor 
receives very little oxygen or nutrients which can lead to 
necrosis [54]. When heat iS applied the surrounding 
healthy vessels dilate while the tumor vessels do not. 
This causes changes in pressure which further decrease the 
amount of blood flowing through the tumor relative to the 
healthy tissue (owen eo Ol. This can lead to _ higher 
temperatures in the tumor relative to the surrounding 


healthy tissue leading to tumor death. 


LeVeen [38] has found that the microvasculature in 
lung tumors is destroyed after a few three-hour 
hyperthermia sessions involving tumor temperatures of 45°C 
Or more. This may be the result of the capillaries 
becoming clogged with erythrocytes. The erumor = tends, | to 
become acidic which can cause erythrocytes to become rigid, 
effectively increasing the viscosity of the blood. This 
effect was discovered and is being utilized by Von Ardenne 
[34]. He increases the glucose level in patients to 


Gecrease their blood pH before applying hyperthermia. 
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METHODS OF PRODUCING HYPERTHERMIA 


Numerous methods have been found to produce whole body 
or local hyperthermia. Each one has its own strengths and 
weaknesses with no particular method best suited to every 
type of tumor. The following paragraphs describe various 
methods of producing hyperthermia and include some of the 


advantages and disadvantages of each technique. 


Whole body hyperthermia has been produced by using 
molten wax baths [34], hot water baths [54], a heated 
enclosure or garment with or without radio frequency (r.f.) 
heating 9§[(54) 11), “injection of “pyrogenic agents [17], 
eee or with heated blood [cited in 54], etc. This type 
of treatment has the potential to destroy unknown or 
widespread metastases if sufficiently high temperatures are 
obtained. A large amount of medical support is required 
because of the danger involved in raising an individual's 
body temperature to the 41.8°C generally used [54, 50]. An 
individual's heart beat increases at the rate of about 10 
beats/°C. This can be fatal especially for those with 
Cardiac problems [50, 54]. Using “Conduction “to. warm a 
person to therapeutic temperatures 1S a long and tedious 
Drocess.)) ne use OlLspyrogens 1s, at, best, (UnSates because 
the response can be difficult to control and unpredictable 


especially when used on a weakened cancer patient [20]. 
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With local hyperthermia high tumor temperatures 
(45=50°C) can be tolerated | =by | patients. Higher 
temperatures mean a shorter treatment time and a more 
erfectivers tumorg, ki li Storm et al] suggest 30 hours of 
Greatmentvat ¢2 °C.) 3 to svshoucs at 45°C, and” a o matter ‘of 


MinuresSsabeoURG to causes thermal deathsofetumors 160). 


Local hyperthermia has been produced by immersion in a 
heated liquid or perfusion of an affected area with heated 
chemotherapeutic drugs, blood, or saline [57, 41, 59]. 
Great care must be taken to avoid permanently damaging the 
heatedveareawl 5/49 138e These techniques have been 
successful in the treatment of melanomas of the extremities 
and of some success in treating bladder carcinomas [57, 41, 
28]. Tumor location and accessibility limit the usefulness 


of these techniques. 


Invasive applicators have been found useful for 
treating “icertain types sof tumors; These devices are 
designed to be inserted into or near a tumor and, using 
microwave or xr.f. fields, produce localized heating [62, 
Dope OC ae Dole Arrays of implantable electrodes have been 
used to localize heating [21]. Implanting metal "seeds" in 
a tumor and heating them with magnetic fields has also been 
tried: for producing! slocalized= *hyperthermiay.[56]% These 
devices are suitable only where the tumor is accessible 
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FUGther,, Implanting, as device in a tumor can produce 
metasteses aS a reSult of mechanical disturbance. Some of 
these techniques are also unsatisfactory because they 


require the use of a general anaesthetic and surgery[36]. 


Ultrasound, r.f. fields, and microwaves can be used to 
non-invaSively produce local or regional hyperthermia. Not 
every tumor is ideally suited to these modes of therapy and 
there can be a problem with producing excessive heating in 


the wrong area. 


Ultrasound is useful because it does not heat fat 
preferentially; temperature meaSurements during treatment 
are simple to make which is not the case with r.f. fields. 
Deep heating iS possible with focussing transducers, 
multiple transducers, or moving transducers [54, 36]. 
Stray leakage of ultrasound into the environment iS not a 
problem. Air is strongly attenuating to ultrasound and the 
mismatch between soft tissues and air is such that 99.9% of 
the energy is reflected at the interface [27, 16]. This 
also means that gas-filled organs such as the lungs or 
intestines cannot be effectively heated with ultrasound. 
The use of a retention enema or the induction of a 
hydrothorax are possible solutions to this problem but may 
not always be satisfactory [36]. There are also large 
reflections produced at bone-tissue interfaces and bone 


absorbs ten times more ultrasound energy than muscle [54]. 
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AS a result, areas blocked by or close to bone cannot be 


treated [16]. 


Radiative microwave applicators operating from 
433-2450 MHz have been used extensively over the last ten 
years. They don't heat fatty tissue preferentially except 
when tissue interfaces create large reflections. The 
applicators tend to be quite compact and the heated area 
can be easily localized. These applicators are usually 
rectangular, OnP Girculanbe waveguides ewith, evor Wit ROUEE 
evelectric Bloadings ore sniages: (26>) 147858, 6425864 | maWith 
surface cooling, waveguide applicators operating at 915 MHz 
have been effective in producing deep heating in human 
thighsSs(s5)eeeThe usew cfs multiple® applicators) “cane 7also 
improve the deep heating ability of microwaves to some 
extent [24, 14, 16]. Unfortunately microwaves have a very 
shallow penetration depth (skin depth 3.7 cm at 433 MHz and 
2 \ecnrateccoONMiziitrelauives tours the ei Skiing depen dy.4 acm 
StmeZien2eMhz and 9le3 cmeat WeMHz) Cor ultrasound energy! so 
they are not generally suited to deep heating[22]. Care 
must also be taken to ensure that resonant effects are not 
BProducedwin’torgans, Cavities, or lambs at the higher 
frequencies [63]. Another disadvantage with microwave 
appl lcatorseuissthatemetale probes | Canes interbere Swath. the 


fields to give false temperature readings [27]. 
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Radiative applicators in the r.f. range have been 
developed. Among these is a ridged waveguide 
dielectrically loaded with distilled water which operates 
ate27. i2)MHz 149, 59). erhissapplicator can be used to heat 
a large, deep area of the thorax particularly when two are 
used in a 'cross-fire' mode [see figure 1.1(b)]. Another 
applicator uses a circle of sixteen transverse 
electromagnetic apertures operating at 50-110 MHz to 
produce deep heating [70]. This particular system utilizes 
water aS a medium to cool superficial tissues and improve 


COuDli ng wl 7 01. 


R.F. energy can _ be coupled inductively to a 
Gieuecieniic7. In this case magnetic fields are utilized to 
induce eddy currents in the dielectric. Often ~a) flat 


pancake coil or a wound coil is used for this purpose [3, 
61]. These applicators do not preferentially heat skin and 
Subcutaneous fat aS much as capacitive applicators do. 
Tissues with high electrical conductivity (je. muscle) tend 
to heat up more than those with low conductivity (je. fat) 


(see page 24). 


The pancake coil produces a toroidal region of maximum 
heat just below the subcutaneous fat layer (Figure 1.2). 
Storm's [61] single winding coil appears to heat an object 
placed inside it uniformly because of thermal conduction. 


Eddy currents circulate around the central axis which means 
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(a) Crossfire capacitive heating 


(b) Crossfire waveguides 
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(c) Rotating electrodes 


Figure 1) ly tModalities tow crossfire, lieating using rst. or 
microwave applicators. 
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Figure 1.2 Heating pattern produced by pancake coil. 
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bherewsise no current produced on the axis [16]..-A plot of 
the H-field inside the coil indicates that the field 
increases drastically aS one approaches the edge of the 
coil [16]. This implies that the strongest heating would 
SGGUn IN the superticialsparts of the bodywif fat-and skin 
hadeedielectric “properties. Similar to “musele. Shaping 
fields is very difficult with inductive heating and metal 
probes can perturb the field giving false temperature 


reagings (275.161. 


Capacitive heating with r.f. has been another popular 
method of inducing hyperthermia. Electric fields are 
produced which in turn’ cause current to “flow in the 
dvelectric. This method is simple and useful because 
electrode size and shape can easily be tailored to fit most 
applications [44, 4]. Multiple applicators produce deep 
heating and the applicators can be cooled to prevent skin 
overheating [37, 65]. Despite the high penetration depth 
of r.f. fields, subcutaneous fat tends to heat up more than 
deep tissues when capacitive heating is used. This can be 
a serious problem when only two stationary electrodes are 
poilized, In addition, the heating pattern can be 
GztEficult tO-predict) and)» metal probes can perturb the rat. 
field giving false temperature readings during irradiation 
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At microwave frequencies phased arrays have been 
suggested as a means of focussing energy deep inside a 
tissue mass without over heating the surface es 253 
Numerous radiators are used and the radiation from each 


phase shifted to produce constructive interference only at 


the tumor site [54]. in practice —this = technique. -15 
unsuitable because a human body is an extremely 
complicated, inhomogeneous, dielectric. As a result, 


attempts to focus energy in this way generally meet with 


failure [14], 


AS previously mentioned, multiple applicators have 
been used to produce deep heating. Generally researchers 
have devised applicators which operate in the 13-110 MHz 
range and have applicators oriented in a cross-fire manner 
(Raicgir emda Cas.b))) (65 > oem ON. These applicators take 
advantage of the large penetration depth of r.f. fields and 
Spare surface tissue, but not the deep tumor site, from 
high radiation levels. As a result, skin and subcutaneous 
fat temperatures can be maintained at acceptable levels 


while deep tumor temperature is increased. 


The cross-fire technique has been used with some 
success by LeVeen to treat lung tumors [38, 39, 65]. These 
tumons are Very difficult. to “heat “by “most —sothner smeans 
because they are deeply embedded in the thorax in an air- 


filled moving organ Surrounded by a bony rib cage. 
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Sugaar and LeVeen have treated lung tumors with their 
cross-fire applicator and then surgically removed the 
tumors for histological examination [65]. These studies 
revealed damaged capillaries in the tumor's stroma leading 
to an immunological response. The immune response was 


Similar to that seen in certain graft rejections. 


LeVeen has had a good success rate, even using a non- 
cross-fire applicator, treating otherwise untreatable lung 
tumors [37, 38]. He found regression in one third of the 
cases treated, relief of pain and improved appetite in over 
80% of the cases, and about 10% of the patients were 
apparently cancer-free for 2-3 years. LeVeen pointed out 
that his cross-fire applicator is vastly superior to the 


applicators used in the above trials [38]. 


An improvement to the cross-fire technique might 
result from using one moving applicator rather than two or 
MOne Switched Stationary applicators [Figure i.1i(c)). The 
moving applicator acts like numerous switched applicators 


BUG WS far ysimpler and cheaper to construcy. 


The use of moving electrodes has already been 
described and patented by LeVeen [39]. He has envisioned a 
system which moves two pancake coil electrodes in a 
semi-circular orbit around the patient. His design suffers 


from using air as the coupling medium between the 


barons 


7 Adhere 4 sa? on 
7  galboal ees 338 ‘dene aif ad 2ai? effiqes Sapenab & var 


eS: ata 
Poatitinn? [ce zipolesgld 7a 


com x 
. Key. stncgees «num «$47 .genegts> sarge [orienta ar 


; py 
trol taggsa Glew" wliedse5 at see Jaks 27 telimes 


= 


202 4 patel 909. ,S476% Feen>54<2 ap 4 668 gad eseViu 

fe. @ Gauss" M6 Jen En! 554 tesiigas sci t+sae 
eg] ta Osid ao ut HO, 82 : ‘ [se .7¢] eae 
sevo ni #ilteqae Garvesqa! Gis fiicy Io lesisy ,bs40039 seRnes 


if Jo roe 
{TPP IAS 


—_ 
= 
ow 
& 
~- 
« 
te 
ve 
~ 
iy 
~~) 
be 
+ 
¢ 
ok, 
LY 
~ 
a 
ay 
> | 


Ki) peed) yhetstd Gee \a@tendewl’ ccioe £) say 


SP. W@ Genglelvne — ote Weaver sd ¢ finelteg Bae bs ies E 
7 7 
ilies ite bie) a sii és vem Aotiy madage 


- 
ea _ mF 


aan Lh a Lite parce thus ‘ a oe ify ‘ bo-t 


. 
- > 


electrodes and skin. This could produce overheating of the 
surface layers. His design also shows electrodes which are 
so” far “apart “that relatively small field gradients are 
present at the patient site. This may have the effect of 
Spreading the fields more evenly over the patient which 
would tend to decrease the advantages of using moving 
electrodes. Further, )hanyslocaliazation ofsfields would: be 
practically impossible and stray radiation would become a 


serious problem. 


For this project a hyperthermia applicator has been 
designed and built which uses moving electrodes as an 
extension of LeVeen's concept. Liisierapplicatorm —utiluzes 
parallel plate electrodes but could use any other kind of 
electromagnetic or ultrasonic radiators. Saline solution 
is used as the coupling medium between the electrodes and 
Skin. Ihis*drastically reduces )rerlectiions (ati athe = skin 
boundary and helps keep superficial tissues cool during 
MoLadiatl. ony, The electrodes are kept within a few 
centimeters of the skin, producing higher field gradients, 
and taking full advantage of the electrode motion. The 
electrode-skin proximity and the large air-saline mismatch 


reduce stray radiation into the environment substantially. 


Immersing a patient in chilled saline would not be 
acceptable as excessive cooling of the entire body would 


likely result. TOqee Coos the irradiated surfaces 
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sufficiently hypothermia could inadvertantly be produced in 
the patient. A doughnut shaped bag filled with chilled 
Saline could be used to encircle the portion of the patient 
containing the tumor. This would require some form of 
thermal and electrical coupling medium between the bag and 
patyvent. The medium would have to be capable of 
penetrating the corified epithelium to be effective. 
Inside the bag electrodes travelling on tracks could be 
used to distribute the rei. over the surface; With this 
type Of “applicator the heating and ©cooling could be 


localized) and readily “controlled. 


Temperature patterns and profiles produced in agar 
cylinders by moving or stationary parallel plate electrodes 
were examined both experimentally and theoretically. In 
each case the surface of the agar was cooled by a saline 
bath. Computer programs HTAVE.0, HTAVE.1, and HTAVE.2 were 
written, and used, to predict moving electrode temperature 
distributions using stationary electrode data [Appendix 1]. 
Stationary electrode temperature distributions were closely 
predicted usSing the computer program MEGAERA developed by 
A. Hiebert Of the University of Alberta Applied 
Electromagnetics Group [30]. MEGAERA simulated, under 
certain circumstances, the temperature distribution inside 
a dielectric exposed to Ieenkee energy from various 
electrodes. With this data one of the HTAVE programs could 


then be used to determine the effect of moving the 
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electrodes. 


Temperature data obtained in this 
demonstrates the Superiority of moving 
Stationary electrodes for producing deeper, 


heating. 
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THEORY FOR R.F. HEATING 


THE QUASI-STATIC APPROXIMATION 


The instantaneous electromagnetic heating in any 
dielectric can be calculated if the electric field 
intensity, E, and dielectric conductivity, 6 , are known 
everywhere inside the Gielectric. The heating 1s 


calculated using the equation, 


P = olE|2 1 


where P is the instantaneouS power per unit volume 
dissipated in the medium. The electric field is described 
exactly by Maxwell's equations subject to appropriate 
boundary conditions. Maxwell's equations are, 


> 
VxE = -yaH 
at 


where oH is. the Magnet cemitdelGgeintensaty jay /aslsee che 
permeabidity, 604s the permittivity, and Jeisecrhe wm current 
density. Unfortunately these equations are very GiherouLt, 


if not impossible, to solve exactly in all but a few simple 
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cases. As a result, approximations are usually made to 


simplify them. 


Consider a medium with 


ieee Bl 3 
WE 
where W is the angular frequency. In this dielectric the 
conduction current is much larger than the displacement 
current. As a result, Maxwell's curl equations can be 


Simoeerred. tot), 


=> 
WxE = -u3l 
at 
> + 
VxH = J 4 


Further, assume that the system has dimensions much smaller 


than a wavelength A in the dielectric where 


Ve 
! 
nn 
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for a uniztorm plane wave [71,50 on pp. 335). In such a 


case the eddy currents produced by the time varying 
magnetic field are insignificant. Hence, Maxwell's curl 


equativens further Simpliry to 171), 
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Tivceeiseethe= —GUaSl=Static= approximation. In etact, these 
equations are true only for static fields [71]. However, 
under the conditions outlined above they can be used to 
approximate time varying fields. Static fields are first 
evaluated then an orale time dependance is assumed to 


calculate heat deposition. 


Equation 6 allows us to define a scalar potential 


d (x2) Such that 


SUDSEIEUcINGs Go wInte 6 Gives, 
Vx(-Vo) = 0 9 


which 1s true because the curl of the gradient of a scalar 


is always zero [33]. Taking the divergence of 7 gives 
V-VxH = Ved 10 


The divergence of the curl of a vector is always equal to 


zero so we have [33] 


Now, 
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which substituted into 11 gives 


V°(cE) = 0 13 
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Substituting 8 into 13 gives 


Ve (O16) a= <0 


Or 
Ve (oVo) = 0 v4 


wheeh “ts a form "of “the current continuity equation. This 
equation is solved numerically by MEGAERA to. evaluate 
fields and temperature Gisioi burr ones inside lossy 


dielectrics. 


Boundary conditions for Maxwell's equations are [33], 
ax(E -E ) = 0 15 
ia 
and 
> > 
ie Bo Er) ) = 6 16 


where E; and §E, are the electric fields in media 1 and 2 
Rots 

respectively, n is the unit vector normal to the surface 

between the two media, wands P.1s theysurface charge on) the 

boundary between 1 and 2% When the Guasi-=static 

approximation is applicable these boundary conditions are 


equivalent to 
= =a 
o E =okE 17 


and 
> = 9 18 
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where O, and o, are the conductivities of media 1 and 2 
— => 

respectively, Ew and ae are the normal components of the 

electric field in media 1 and 2 respectively, and @ and 

?, are the potentials in media 1 and 2 respectively [30]. 

Equations 17 and 18 are used by MEGAERA when solving the 


Cur rentmconramnuiey equacion 41.350 )> 
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HEATING WITH COILS AND PLATE ELECTRODES 


Tissues with high electrical conductivity (je. muscle) 
tend to heat up more than those with low conductivity (ie. 
fat) when inductive coupling is used. This occurs because 
the applicators produce E fields which are usually parallel 
to the fat-muscle interface. Boundary conditions at the 


interface are such that, 


E.= 19 


where Bes is the tangential field in the fat and a as the 
tangential field in the muscle [figure SAE The 
conductivity of the fat is lower than that in the muscle so 
the amount of power deposited in the fat is lower than that 


in the muscle. This can be expressed as, 


= Esso 20 


where P. is the power deposited in the fat, P,, is the power 
deposited in the muscle, O, is the fat conductivity, and 


Semis the Musclenconductivi ty. srromeediation | oewe have 
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PANCAKE COIL APPLICATOR 
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Pigure, 2.1 E - fields inductively coupled into tissue byecoul 
type hyperthermia applicator. 
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The rationale for the high heat production by 
Capacitive applicators in relatively low loss fat has been 
Simply illustrated, as follows, by Christensen and Durney 


[16]: 


Consider two long capacitor plates with fat and muscle 
dielectric layers between them (Figure 22 From 


Maxwell's equations we have 


jut, ; — : 
where an e time dependence 1S assumed for E, and Gis the 


aielectric conductivity which contains all loss terms. 


Taking the divergence of both sides gives 
Vs (Gat jwe)E = 0 


Using Gauss' theorem and the usual pillbox argument ,across 


thes boundary in Figure 2.2 5[5ip. 95), 10 ‘can be shown: that 


E E 22 
Secnt a oan 
Bue iS the normal r.meseelectricl field in the fa tebe is 


Eheunonmal rms. electric field in the musele <f s is othe 
permittivity in the fat and €, is the permittivity in the 
muscle. In the fat the average power absorbed per unit 


volume Pe is 


Pulsioe bone 23 
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CAPACITOR PLATE 
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CAPACITOR PLATE 


Ricures2 <2 — - fields capacitively coupled into tissue by 
parallel plates. 
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where GO, 1S the conductivity in the fat. Similarly, in the 


muscle P. is 


(4 24 
where 9, is the conductivity in the muscle. The ratio 


2 SS 25 
=>: 


describes the amount of power absorbed in a unit volume of 
fat relative to a unit volume of muscle. From equation 22 


we have 
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At a frequency of 27.12 MHz we have for fat and muscle 
Em =113-3405.82, €¢ =20.0-57.23, G% =1.09rI0S/m, and 0,,=0.612 


S /m [tore Thus, 
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However, 
JE gl? | 113-9405.82|? 
ee = 502 8) 
Eee | Sw 2020 =47,223 = 
mm 
Therefore, 
°% 
—— = (0.02) (92.37) = 7.85 
P 
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Or, there is more heat in the fat than in the muscle. 
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MEGAERA: THE FINITE DIFFERENCE APPROXIMATION 


The following section, content and notation are drawn 


from the work of Hiebert [30]. 


MEGAERA uses a finite difference approximation to 


solve the current continuity equation 


Ve(oV¢) = 0 Bi 


In two dimensions this equation is: 


@ G0 
2-28) + = 22) = 0 32 
ox sox dy dy 


For the first term in eqn. 32 the finite difference 


approximation 2s% 


where =O, and | ¢) ateutherconductivity and potential vat the 

mtn : ; stn 

zt POCaeLOn ..anG hj is the distance between the 1 and 
th , ee 

Piel location 8 ihigure e243) ’evAGsimilartequation 


exists for the second term in eqn. 32. 
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The area of interest is divided up into numerous 
rectangular blocks with centres defined DY i, 
coordinates. Each block is assumed to be completely 
homogeneous with regard to electrical and thermal 
properties. If one applies the two dimensional finite 
difilerence equation to grid block number (i,4) and the four 


neighboring grid blocks [figure 2.3] an equation of the 


form 
Re 2G, ph tS)) ab: eee at, ee = 
i,g°i,jtl i,j°itl, 3 i,g°4,3 Pe 11rd eae a a 
results, where 
pissy c ae tinted ga jagan? 
= ate 
Sig = Spar, 74g h Oty)? 
A. , = -R, .-S, .- =F 33 
264) 1,J 1,J) 1,J) 1,J 
= + 
C. (O54 tg 5)/ ty, (hy ee 
2 + 
Fi (Fp pty /tk;_1 &; k5_4)} 


and k; is the distance between the (i,j) and (i,j+1) 
positions. These five equations, 35, in five unknowns are 


solved to find the coefficients in equation 34. 


An equation similar to eqn. 34 1s found for each) (grid 
block in the domain. The resulting system of equations 
forms a matrix which MEGAERA solves with the aid of an 


iterative relaxation method known as the Alternating 
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Direction Implicit Procedure. 


At the boundary between an insulator in grid (i-1,3) 
and gyamedielectric sine grid. blocky c(i), 5) the boundary 


condition, 


ay 
= 0 36 


where n is the normal, is satisfied by setting 


On Oe 37 


Pie chic block (i=l) j) 1s an electrode then Oe. is known 
and a modification is made to its coefficient to allow for 
the assumed infinite conductivity. When grid block (i-1,3) 
nas a ditferent conductivity «than block 41,5), the boundary 


eonditions 


1 M1 22 


and 


aresutibized=townodiufy thescoefficients. 


As previously noted, the average heating rate produced 


by a sinusoidally varying electric field is given by 


P = of B|? 40 
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When evaluating the temperature distribution inside a 
dielectric the only heat transfer mechanism taken into 
account by MEGAERA is conduction. While convection and 
radiation play a part when heating an immersed agar 
cylinder their effect is small. Infra-red radiation does 
not penetrate far in water and convection does not take 
place in the agar gel at the temperatures used. The water 
bath is sufficiently large that it can be approximated by 
an infinite heat. sink. Water jets are used to direct a 
flow of saline over the agar surface so convective effects 


between the agar and and water bath can be ignored. 


The equation used by MEGAERA to model heat production 
and flow is 
MOT _ Ve (kVT) + o|E|? 41 
dt 
wiere M 1s the volumetric) heat capacityein J7°Cm°) k a6) the 
thermal  CONGUuctiVity s1n W/meCc, Tus the temperature: yn C, 
aid tb eis stvne= im “seconds. The finite difference 
approximation for the spatial portion of eqn. 41 1s similar 
fovthat for the Current Continuity equation. Por | example, 
the term 
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Boundary conditions are handled in a manner analogous to 
the electrical case. The time derivative in 41 is 


evaluated using the forward difference approximation, 


MOT . M, Hee ee 
aye antec Soa 43 
At 


where n refers to the number of the timestep and Atis the 
Size of the timestep. The heating rate term in eqn. 41 is 
calculated using the current continuity equation for the 


finite difference equation. 


The thermal finite difference equation 1s solved by 
using the.» Alternating sUurecti on MIMmplncit eMethod-=s This 
method is used for parabolic partial differential equations 
while the previously mentioned Alternating Direction 


implicit Procedure Ys used ftcmsolve ellipive equations. 


Using the above methods and approximations Megaera is 
ablestousimulate 1.f. Neatang. The «power Of the: program ais 
its ability to simulate the heating produced by any shape 
or combination of electrodes. The accuracy of this program 
is demonstrated in following sections and used extensively 


to model different electrode configurations. 
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METHODS AND MATERIALS 


VOLUMETRIC HEAT CAPACITY MEASUREMENT OF AGAR 


The temperature of a (2.6+0.05) cm diameter’ by 
(3,640.05) cm long agar cylinder was measured then placed 
ina thermos: containing § (230.023.0) gm ‘of "heated water 
(Figure 3.1). The thermos was capped with styrofoam which 
was securely taped in place. A wood handled stirring 
paddle and a thermocouple extended through the styrofoam 
Cape anto “the heated, water. A Fluke 2190A Digital 
Thermometer was connected to the copper-constantan 
thermocouple. Water bath temperature was monitored over a 
period of time to determine the amount of energy absorbed 
by the agar and the amount lost through the calorimeter. 
With the volume of agar and amount of energy it absorbed 


known the volumetric heat capacity Cyscould be calculated. 
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STYROFOAM 


FLUKE 2190A 
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HEATED WATER 
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Figure 3.1 Apparatus used for heat capacity measurements. 
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Thermal Conductivity Measurement 


To calculate thes temperature distribution in rt. 
heated agar the thermal conductivity k had to be measured. 
This was done using a measurement chamber belonging to the 


University of Alberta Applied Electromagnetics Group. 


TWoseloeetOns peCm maOvaMerer Cylindrical. sections | Of 
agar dielectric (0.9% NaCl, reagent grade, and 3% agar by 
weight dissolved in distilled water) were Sandwiched 
between three brass disks (Figure 3.2). The upper agar 
SlabewassvG|.iot0. 03) echemonick © and the lower Slab 
ie 55+ 00s) mcme thick The upper brass disk contained an 
electric heater coil and a copper-constantan thermocouple. 
Between the two agar slabs was a thin brass disk containing 
a thermocouple. The lower brass disk contained a 
water-flow heat sink and a thermocouple. Contact surfaces 
between the brass disks and agar were lightly coated with 
Silicon Thermal Compound #120-8 (Wakefield Engineering 
Company) to ensure good heat transfer. Plastic™ “trim | was 
wrapped around the sides of the agar slabs to retard 
evaporation. The agar and brass disks were insulated by 10 
CMoemon styrotoOam, tOpsand bottom, and a minimum of 13 .cm) of 
styrofoam on the sides. For added support and insulation 
the whole chamber wasS contained in a tightly fitting 


plywood box. 
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Figure 3,2 Chamber used for thermal conductivity measurements. 


Power waS Supplied to the upper brass disk via an 
HP6439B D.C. Power Supply (Figure 3.3). A Keithley 168 
Auto-ranging DMM and a Universal AVO Meter (10 amp. scale) 
were used to measure voltage and current respectively. A 
Fluke 2176A Digital Thermometer was used to measure 


temperatures. 


The middle brass disk allowed temperature measurements 
to be made between the heat source and heat sink. In this 
way three thermal conductivity measurements could be made 
at the same time. MeaSurements of k with the thin agar 
slab, the thick slab, and both slabs together compensated 


for convection along the sides of the agar. 


The Fourier heat conduction equation [27] 
O=kTACE =i) /x 


where, Q-input power (Watts) 
A-surface area Of agar cylinder top (m7) 
x-thickness of the agar slab (m) 
T,-temperature of the heated agar face (°C) 
T,-temperature of the cooled agar face (°C) 
k-thermal conductivity (Watts/m°C), 


was used to evaluate the thermal conductivity of the agar. 
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Figure 3.3 Apparatus used for thermal conductivity measurements. 


Dielectric Measurements 


The dielectric properties of the saline bath and agar 
cylinder were determined using a method described by 
Bottomley [7}). A nylon disk, 6 cm in diameter by 1° cm 
thick, had a 1cm diameter hole drilled through its centre 
(Figure 3.4). The hole was countersunk on either end by 
ecncentric, 2 cm dyameter, holes 0.35 cm deep, Two small 
holes, 1 mm in diameter, were drilled in opposing sides of 
the nylon disk, through to the central chamber. These 
holes were fitted with 19 gauge hypodermic needles and _ 10 
Comms yringes: The meaSuring cell was oriented so that the 
Syringes were at the top and bottom of the disk. The lower 
Syringe initially held the fluid to be measured and the top 
Syringe acted as an overflow reservoir and drain for air 
bubbles. Non=ligquid samples) had) to be carefully cut to 


Suze stor tne cell, 


A silver plated copper disk, 2 cm in diameter, was 
placed in e€ach countersunk hole and held ain place “by 
Styrorcame Locks mand mplastic —tle-stlaps. The area of 
contact between the nylon and electrodes was made 
waterproof with Dow Corning High Vacuum Grease. Short 
wires (5 mm long) were soldered to the electrodes and 
connected through an HP00601A component mounting adapter to 


an HP4815A Vector Impedance Meter (Figure 3.5). Frequency 
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Figure 3.4 Dielectric measurement cell. 
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Figure 3.5 Apparatus used for dielectric measurements. 
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was monitored with a Fluke 1900A Multi-Counter. The sample 


temperature was monitored with a copper-constantan 


thermocouple, in contact with the ground electrode, 


Fluke 2100A Digital Thermometer. 


a 


The conductivity and relative dielectric constant of a 


Sample in the cell can be shown to be [Appendix #4] 


rae zcosd 
OEE ee ee 
CZ Sea) Le) 


and 
WL -zsind 1 
cy A ey ti zsind+tw-L ae Sees) srl 
P Pp O 
where, 6-the sample conductivity 


z-the magnitude of the measured impedance 
$-the phase angle of the vector impedance 
2 -the cell constant 

W-the angular frequency 

L,-the probe inductance 

Co-the empty cell capacitance 


C,-the probe capacitance 


€ ~the sample relative dielectric constant 


€,-the permiuutivavy sot free spacer 


The cell was calibrated by measuring the impedance 


of 


low conductivity liquids, with known dielectric properties. 
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A frequency of less than 5 MHz was used so that L (about 10 


nH) could be neglected. In this case the simplified 


expressions 
cos 
il RZ 
and 
~ sing il 
ig WZ re COT ae R 


fe) 


could be used. A graph of €, vs. -(sind)/uz was plotted 
Givingva straignt Wine’ withe a Slope of (ee) eeats a 
V=Intercept (Of (Go nn) cer (Figure 3.6 and Tables 3.1 and 
3.2). Values of (€, £)7'=6.11x10'?F and (Co+C )=4.73 pF 


were obtained. 
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FIGURE 3.6 


Methyl alcohol 


© Ethanol(95%) 


Isopropyl alcohol 


Slope=l__=6.108x10!29/s 
(ed 
© - 
£=1.849x10 2m is the 
cell constant. 
yointercept=-{C 4¢ ) 
ef 
=-28. 866 
+C_)=4, 
ae oe 4.726 pF 


=the capacitance 
of the HP probe plus 
the empty cell. 


10 


47 


plas Tew cae 04 a ied 


do vidinas % 
% 
\ : * 
ieee Ge ave x 
wT 
Dep! = 
a i ai a 
- (nin a me 
es @FA Wee § 
ee a 7 
——— 7 oa 
@ 
ea (adie Mt? 
eu Vivet9 
7 


TABLE 3.1 Calibration Fluids For Dielectric Measuring Cell 


DIELECTRIC FREQ. (MHz) e _[30]* ne (x10 12) 
(s/n) 


METHYL ALCOHOL 


ETHANOL (952) 


ISOPROPYL ALCOHOL 


ca DMS Cs 
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Temperature Distributions 


An agar cylinder (0.9% NaCl, reagent grade, and 3% 
agar, by weight dissolved in distilled water), 10 cm in 
diameter by 3 cm thick, waS immersed in 300 litres of 
Saline solution (0.9% NaCl, by weight, dissolved in tap 
water) contained ina large polyethylene tank. The saline 
and agar were found to have similar dielectric properties. 
The conductivity of the saline at 3.75 MHz was measured to 
DemrOn=(2.0)20205)5s/MeCremagaiep the picondvctivity Gwas 

Orne eo oL0 02) S/m: From published literature the 
Felatave, permittivities Of theragar)) and» saline, “at 3.75 
MHz, were estimated to be €, = €, = 8042 [31, 32]. The 
Similarity in electrical properties permitted good coupling 
of the r.f. to the agar and decreased the sensitivity of 
fof absorption to any inregulamstres in) the “agar. shape. 
Saline prevented salt leaching from the agar, kept the 
Gielectric homogeneous, and acted as a heat sink to cool 
the agar surface. The cylinder was thermally insulated on 
its top and bottom faces by styrofoam caps which reduced 
heat loss except in the radial direction. The caps were 
wrapped in plastic film and sealed with Dow Corning 
Silicone Sealer to keep water absorption from compromising 
their insulating ability. A plexiglass clamp held the agar 
cylinder and styrofoam caps securely below the saline 


SUibacem rl Ores.) ). 
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Figure 3.7 Apparatus used for heating agar cylinder. 


50 


i 
» y 
a aorta, 


| : hea UT Sy i) 


Seay 


The irradiating electrodes were each 1 cm wide by 15 
cm high and made from flat silver-plated copper. Tall 
electrodes produced an E-field in the 3 cm high agar 
cylinder. The E-field did not vary greatly in the axial 
direction (Appendix #3). This permitted two dimensional 
Simulation of the heating to be made. The electrodes were 
Supported by nylon posts and could be made to revolve 
around the agar in an oscillatory manner, or held 
Stationary. Amplitude and speed of oscillation as well as 
electrode distance from the agar surface could be varied 


(Figure 3.8). 


R.F. power was provided by an HP615A Test Oscillator 
and amplified by an Electronic Navigation Industries Inc. 
PNAS Oar ry. Power Amplifier. The oscillator was 
powered through a GraLab 600 Timer. A Decca KW-109 
Supermatch was used to match the amplifier to the load 
(Figure 3.9). 'Low' frequency matching was optimized by 
Bie eekKWo109meana. the selectrode configurations used, at 
(3.75+0.02%) MHz. Frequency was monitored with a Fluke 


1900, Multi-Counter. 


Temperatures were meaSured at nine points in the agar 
and one in the water bath. Copper-constantan thermocouples 
(1 mm diameter), a Fluke 2190A Digital Thermometer and a 
Fluke Y2100 Thermocouple Selector were used to make 


temperature measurements. The thermocouples were oriented 
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Figure 3.8 Apparatus used to move electrodes around heated 
medium. 
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Figure 3.9 Apparatus used to maintain saline bath temperature 
and heat agar cylinder. 
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perpendicular to the E-field to minimize current induction 
and artifactual heating. They were implanted to a depth 
mid-way down the length of the agar cylinder and along two 
radial lines perpendicular to each other. One line was set 


parallel to the stationary electrode E-field CRA Cuners. 9) 


The water bath and agar were maintained at a 
temperature of (22.1+0.1)°C by a Brinkman IC-2 water bath 
pump/heater. The Brinkman was controlled by the Fluke 
2190A through a relay, forming a feedback loop. A water 
bath temperature just above room temperature was chosen to 
avoid the need for a heater/refrigerator unit. Water jets 
located beside each electrode directed the flow of saline 


over the surface of the agar to aid cooling. 


With thermocouples in place and the digital 
thermometer switched off, the agar was irradiated for 10 
minutes at an output power of 200410 Watts. Wuls(ene oe nea 
amplifier was then turned off and all ten temperature sites 


monitored over a period of about 6 minutes. 


A curve Piece computer program was used to 
extrapolate the cooling curve data for each ‘thermocouple 
back to the time the r.f. was terminated [23]. The program 
Gave various Statistics) to) indicate: the reliability sof veach 
Extrapolabedmapoint ma(hIgUrem~s..0—sand | Table sesa2)., The 


curve-fitting program was used to avoid 
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COOLING CURVE WITH 99% CONFIDENCE BANDS 
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POCURE 9. LO 


TABLE 3.2 Fitted and Measured Cooling Curve Points 


Pi Pepe. RESIDUAL MEASURED Y 


24.550949 0.049049184 
24.315445 -0.015454937 


55 


of c.3ee ) iis os . ra, 


dase ora 
ion " ee Gael ee _ 


of v.84 


any bias in plotting the cooling curves. This procedure 
was repeated for several different orientations between the 
electrodes and thermocouples. Then the agar cylinder was 
cut in half to determine more accurately the position of 
each thermocouple. These data were then used to plot the 
temperature distribution inside the agar cylinder after 10 


minutes of exposure to the r.f. 


Data points obtained as above resulted ina fairly 
Goanrse —emperalurerdistribuuion (Figure © 3.4 lps. Plots of 
temperature vs. angle were made for various radii from the 
Centre of the “agar phantom = (figures 3.212)% From these 
curves temperatures at points which were not measured could 
be estimated. In this way a denser field of temperature 


Gata could be obtained for improved plot production. 


To predict temperature distributions resulting from 
moving electrodes, computer programs were written which 
used static-electrode temperature distributions as data 
(Appendix HAD. These programs calculated a weighted 
average temperature over a specific radius and arc length 
ipsaoe. the agar. Radia of —0, 1;275,4, and-4-05  cmewere 
used. The arc length corresponded to the orbit of the 
electrodes and the calculation considered the variation in 
temperature along its length. FOr the measurements 


considered here the arc length was +40° or 80° total. 
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The weighting functions used to calculate the average 
temperatures were based on the electrode motion. This 
resulted from the fact that the temperature increase at a 
given point was inversely proportional to the speed of the 
electrode past that point. In other words the temperature 
increase varied proportionally with exposure time to a 


Particular k-field: 


Assuming that the electrode speed varied sinusoidally 
from 0 m/s at the beginning of its swing, to some maximum 
at the centre of the swing, and back to 0 m/s at the end of 


its swing the following weighting function would result: 


T=07 04707, 7050501014477, )+0,0614 


(T3700), OS41( Ts 1s) +0. 2709T, trae 


where, T-the averaged temperature, in °C, which would 
occur at. the centre of an 80° anc 
T o-the temperature, in °C, measured at the centre 
OL tne warc.. 
D.—the temperatures, im °C, measuredvat = (1 0xn)° 


respectively, from the centre of the arc. 


The constants in equation 44 are inversely 
proportional. to. thes velocity ot che electrodessover = that 


pointeim the arc. | The sumsot the constants over sone entire 
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are as normalized to “unity. By applying the weighting 
function to the static electrode temperature data a plot of 
temperatures (which would result from electrodes with a 


Sinusoidally varying speed) was obtained. 


Water resistance and play in the electrode moving 
assembly prevented a true sinusoidal variation in electrode 
Speeqc®) from soccurnring,. The electrode speed variation 
elosely resembled a flattened sinusoid or half circle. For 


the circular speed variation the weighting function was 


T=O0700)1 Tee 00832 (Cl +T, +08 0905 


(Tey ls Olt OC oMmist Ie) tO (ONT Dae) 


The flattened sinusoid was approximated by averaging the 
Circular and sinusoidal functions. This weighting function 


was 


T=0.0644 7570 ..0666( 7,47. )+0.0760 


Cle eis) Orem Os4 (Cle wren len? OO) (ee ean) 


Temperature distributions calculated using the 
weighting functions were compared with measured data. The 
Gata wes measured in the agar along 0°, 45°, 90°, and —135° 


radial lines (Figure 3.13). 
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Ereure 37135 |Top view of Jagan cylinder showing radial. lines along 
which temperatures were measured. 
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Great care had to be taken to ensure that the agar 
phantom did not move. during a Set of measurements. 
Further, the phantom had to be carefully positioned so that 
it was concentric with the electrode axis of rotation. ite 
the phantom moved off centre the stationary electrode 
heating data would not be useful. As the electrodes moved 
their distance from the phantom would change. This would 
cause a change in the E-field with electrode position which 
would invalidate the technique used to. produce static 
electrode temperature maps. As a result, only temperature 
patterns obtained from electrodes moving along a path 


concentric with the phantom were used. 


Non-concentric electrode motion could have been 
considered if provision was made for rotating the phantom. 
Care would have to be taken to ensure that the 
thermocouples did not move inside the agar. Complete 
temperature maps of the inside of the agar would have to be 
made for many electrode positions. These maps could. then 
provide the data necessary to predict the moving electrode 
temperature distribution. Computer simulation of the above 
case would require mumerous Megaera runs to simulate the 
various electrode positions. This entire amount of data 
would then be averaged together and weighted according to 
the electrode speed to produce a moving electrode 


temperature map. 
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The non-concentric electrode motion could prove useful 
with regard to moving and shaping an internal hot spot in 
agar. It was not examined in this work primarily because 
of the increased cost in computer time that would have 
resulted. In addition, the geometry of the problem 
indicates that concentric electrode movement would result 
in a higher centre temperature to peak temperature ratio 
than non-concentric motion. Further, the inhomogeneity of 
humans would be such that accurately shaping and steering 
an r.f. induced hot spot would be nearly impossible. So 
the concentric motion case would be of greatest interest 
for producing deep heating in a clinical situation, without 


overheating superficial tissues. 
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Live Animal Measurements 


A 3.5 kg male cat waS anaesthetized (Nembutol) and 
Shaved around the lower half of its torso. A midline 
incision was made in the ventral surface and the right 
kidney exposed by dissection. The right renal artery and 
veim were tied off to stop blood flow to the kidney. In 
this way the kidney roughly modelled a solid tumor. A stab 
wound was made just below the animal's right shoulder 
blade. Copper sheathed copper-constantan thermocouples 
were passed through the wound to _ the kidney site. 
Thermocouples were embedded in tissue along a line 
intersecting the centre axis of the cat and the right 
kidney. Thermocouples were also placed along a line 
perpendicular to this. AS much as possible, the 
thermocouples were kept perpendicular to the E-field. The 
incision was closed using 3-0 silk and waterproofed with 
vaseline. A cylindrical cloth harness held the animal 
upright with, its head above the saline.” The harness “held 
the cat between the rotating electrodes and kept the 
animal's torso cylindrical. The electrodes were’ each 2 cm 
wide by 4 cm high and made from silver plated copper. The 
animal wasS approximately 12 cm in diameter and the 
electrodes about 4° em@irom its skin. To avoid producing 
hypothermia in the cat the saline bath was maintained at 


Tim lod Get The animal was irradiated for 20 minutes at an 
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Output power of (250410) watts. Electrodes orbited around 
the cat in an oscillatory manner through an are of 80 


Temperatures were meaSured as with the agar. 
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DATA AND RESULTS 


MEASUREMENT OF VOLUMETRIC HEAT CAPACITY OF AGAR 


Data obtained from volumetric heat capacity 
measurements are shown in Table 4.1. The water was 
inca eat (66.320 ..1)@C and “the agar “cylinder. at 
(20.1+0.1)° C. When the agar was placed in the calorimeter 
the water temperature began decreasSing exponentially 
[Figure 4.1]. After 1000 seconds the temperature began 
decreasing linearly. The exponential decrease was the 
result of energy absorption by the agar. Energy leaking 
from sche) calorimeter!’ produced the linear “portion. “By 
extrapolating the straight portion of the curve to t=0 
seconds the loss of heat through the calorimeter could be 
accounted | for. Figure 4.1 indicates that the water 
decreased (3.5+0.2)°C as a reSult of the agar, conversely 
the agar increased (42.7+0.2)°C as a result of the water. 


The total amount of energy absorbed by the agar was 


OP ar =(3.520.2) (Cx 230.073 50) qm=(805t33.0) cal 
=(805t33.8)calx4d. 164 J/cal=(3368 15141 04)3) 


Heat capacity is the amount of energy (J) absorbed by a 
Material, divided by the change in temperature (°C) of the 


material. So, the heat capacity for the agar was 


Ge (2sbeul le toy em e0. 2). C= (er Oren) Jac 
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TABLE 4.1 Water temperature in calorimeter versus time. 


TIME (seconds) TEMPERATURE (40.1 °C.) 
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The volumetric heat Capacity, used by Megaera, is the heat 


capacity of the material per unit volume (m?). For the 


agar sample used, this gave a volumetric heat Capacity of 


C=(78.943.7)5/°Cx(5.2+.03)x10 ‘nt? 


=(4°4+0.4)%10°9/°C m: 
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THERMAL CONDUCTIVITY MEASUREMENT FOR AGAR 


From the Fourier heat conduction equation 
k2=(0 MO Ay Te Lev x) AAR a Tol) 


Pomsthe il 1520.02 )emiehich slab Cf.agan,. using the data in 


Table 4.2, 
keGOn52+0 303) W/m G 


Table 4.3 shows values obtained for different thicknesses 


Ob agar. 
Inaccuracy in the k measurement increased with 
increasing slab thickness. This occurred because energy 


was more readily lost through air convection around a thick 
Slice of agar and imperfect insulation. Unfortunately the 
precision with which the thermal conductivity could be 
calculated was inversely proportional to the agar slab 
thickness. This wasva result of uncertainty in ) the slab 
thickness. To arrive at a reasonable approximation for the 
thermal conductivity of the agar, regression analysis was 
WMSeameEO wind tne thermaleconductivity fOr am 0s Oscm um thick 


Slab (Figure 4.3). 


The straight line graph in figure 4.3 was determined 
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TABLE 4.2 Temperature measured at three sites in agar cylinder 
for given voltage and current . 


Voltage = (3.505+0.005) V. 
Current = CIE 720-01): A. 


* See Figure 4.2. 


TABLE 4.3 Thermal conductivity, k, measured for various 


thicknesses of agar. 
T)-T> (CC) | k Gi/m-C) 


AGAR THICKNESS (cm) 


1.520508 Peace OL OD 220505 
Sess Uae) 2 Oe Oem o2 0015 
AoOeO.0s 26000 1 OS Sea Git 
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Figure 4.2 Arrangement of materials and apparatus inside the 


thermal conductivity measurement chamber. 
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where n=3 and, for a line of the form y=at+bx, the terms 


w,Y,x° and xy given in) Table 4.4, From the above, 
a=0.51 

and 
b=0.0098 


Taking the worst case limits, as shown in figure 4.3, the 


Ehermar Conductivity for agar was found: to be, 


kK=(0.5120 04) w/m°C. 
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TABLE 4.4 Values used to plot straight line graph of agar thermal 
conductivity versus agar thickness. 


es) Oe52 a2 0.60 
Sess) Oeo> 122,60 gs) 
4.70 Oro 22709 
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THERMAL CONDUCTIVITY (W/m°C) 


CALCULATED THERMAL CONDUCTIVITY FOR VARIOUS AGAR SLAB THICKNESSES 


AGAR THICKNESS (cm) 


FIGURE 4.3 
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ELECTRICAL CONDUCTIVITY MEASUREMENTS 


Saline solutions of known molarities and dielectric 
properties were used to determine the accuracy of the 
dielectric meaSuring cell. The cell proved accurate for 
measuring the conductivity but not the permittivity of 
highly conducting solutions [Table 4.5]. This was most 
pronounced at lower frequencies and in solutions with 
higher molarities. (hejsource fot this error, -lay “inthe 
inability of the HP 4815A to accurately measure very small 
phase angles. This led to a large amount of uncertainty, 
or lack of precision, in the permittivity measurements. 
The relation used to calculate the permittivity was far 
more sensitive to changes in a small phase angle than was 


the equation for tconductivity. i.e. 


o 45 


46 


The reason for this can be seen when the sine and cosine of 


two small angles, like 1° and 0.5°, are compared.i-e. 


Samii. Or v=0 01745 
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Sim. 5°) =0. 00873 
cos(1.0°)=0.99985 
cos(0.5°)=0.99996 


A 0.5° change in angle produced a 50% change in the value 
of the sine function. The same change in angle produced a 
0.01% change in the cosine function. Clearly, equation 46, 
with its sine dependance on bd, should be more sensitive to 


Variations in small angles than equation 45. 


The water bath, used in the r.f. heating experiment, 
hac mam MOlalsiiiy On suaDoUren O01 5 and the SUgnuiwocant 
measurements made at 3.75 MHz. This combination produced a 
phase angle which was unacceptably small for accurate 
permittivity meaSurements. As a result, extrapolations of 
published values were used to eStimate the relative 
permittivity of the water bath, and agar, as 81+2 [31,32]. 


Conductivity measurements were quite acceptable. 


The conductivities measured for the agar and water 


bath proved to be very consistent giving average values of 


al oad 
Ons (22 0120-03 )s) =m 


and 


oO (1.99+0.03)2 !m } 
agar 


These two values would be expected to be very similar as 


Od GyT3L0m2 6 
7 ae whan het esleng Ciena mie 


pianel regen Qeiidedt’ »3.> gre nt) Sean inte wee 
wrenbripgis «6 ald, «Bas 21 ¢ wate, Ta. ephsetee 0 (be 
gf beeitorg yebjanbieas a22T 144 GF .G iw aba | 
/qhesinds oa bese witwalass vey, Bev” woRtye signs 
‘ ie stialialoniziap @hite? 5 @ qucsmbsirenn (ystebseb 
idaler aif? \“éeppaea cf Mh cy etwisv ae 
AR 7h) c27e oe | see Site tect vee ei? Se nisl 


elhagendeg Vitlup Oo 1c Si eIegse qieetanunes 


san¢pehz>ubeang .66F 


’ 


Zateau 8F4 ipa Wt aed | i“) © 
te AGUMAY BEA Xe eagle in’ itest foo GaIey. AS oa Gavesq haed 


ri deiy 


eel eg-inie sh) -- 
sane ® 


Me FAT WOER YY, Els Dervaus? oF DLiee outtay ae? 


& 


nn | a 
“ss ae : 2 


PONSMCOULAS CaSil yee ciEruscmuintO wssOr OUts Of the agar 
SbIanging, it) into sequilibrium with the “bath iio2).) The 
measured conductivities were high for normal saline 
solutions. This was expected as tap water was used in the 
bath. Further, the experiments were performed during early 
Spring which, in Edmonton, means the water had a higher 
impurity content than uSual. As long as the conductivities 
of the bath and agar were meaSured the problems posed by 
uSing tap water were deemed acceptable. Tis was 
especially true in view of the difficulty in obtaining and 
maintaining large quantities of distilled water. AS an 
added bonus colonies of bacteria and mold were unwilling to 
grow on agar immersed in the tap water. With distilled 
water in the bath black patches would become visible on the 


agar after a few days. 


For very poor conductors the large phase angle would 
lead to an unreliable conductivity measurement and an 
accurate permittivity measurement. Dielectric measurements 
on materials, and at frequencies, which gave a phase angle 
around 45° would produce the most accurate conductivity and 
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From the above data the loss tangent for the agar was 
tanS= 11745 . 


The wavelength of 3.7/5 Mizar ot.sin the agar Was 
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With tan 6 ey cn Gaset hewsvagar Cylinder much less than a 
wavelength across, the quasi-static approximation could be 


used, 
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TEMPERATURE MEASUREMENTS IN THE AGAR 


Ones Source SOlesuncertainty with the temperature 
measunemenuS §WaSmethat ms resulting [romecurve siuttingeand 
extrapolation. The statistical error associated with 
fitting curves to the cooling data and extrapolating back 


to the time the r.f. was terminated was evaluated using 


[S51 


2_y» 42 eats 
Lang, Les 
Our Merarogeg) (log ata) 


where, 
y-the measured temperature, 
y-the temperature predicted by curve fitting, 
x-the time of each measurement, 
x-the average of the times the temperatures were 
measured, 
n-the number of data points obtained for each 
curve. 
TypICouLy es othits Bernore Never exceeded 07 oC for anys daca 
point (Figure 4.4): Drift in the power source (45%), 
VAarLatiOons © an agar radius (23%), and error imherent ingche 
use of copper-constantan thermocouples (1+0.1°C) were three 
other majour sources of error. Maximum drift in 
temperature data taken from the centre of the agar phantom 


amounted to no more than +4%. This was more indicative of 
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the reliability and accuracy of the measurements than the 


sum of all previously mentioned uncertainties. 


Drift in the digital thermometer was insignificant 
over the course of the experiment. The resolution of the 


Fluke 21490A was limited to +0.05°C. 


The water bath temperature was controlled to within 
+0.1°C for each experiment. The volume of saline was great 
enough that no temperature increase due to r.f. absorption 


was detected. 


Uncertainty in thermocouple position was rendered 
instoni facantubystateraily slicing =the» agar scylinder =to 
find the thermocouple holes. This could not be done with 
ENewecaes. Handling of the animals would likely have 
shifted the thermocouple positions a considerable amount. 
This uncertainty was estimated to be about +0.5) cm for most 


points. 


In eaddition to the above uncertainties, the 
measurements made in live cats suffered from another 
problem. Thermocouples, in), the “cat "were not always 
perpendicular to _ the E-field, hence an unknown amount of 


error could have been introduced. 
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MEASURED TEMPERATURE DISTRIBUTIONS 


The temperature distribution produced by stationary 
electrodes, each (1295+0.05)cm from the agar surface,’ is 
shown in Figure 4.5. Figure 4.6 indicates the relative 
tocation of the electrodes to produce this” distribution 
plus x and y axes for reference purposes. As expected 
there are two symmetrically placed temperature peaks, one 
near each electrode. The peaks occur about 1 cm deep 
inside the agar because of the cool saline bath. If the 
Saline was replaced by a thermal insulator the temperature 
peak would be much larger, occur at the agar surface, and 
decay exponentially toward the centre. Clearly, the 
advantage of using a cool bath or bolus iS apparent for 


Gooling superficial Gielectric or tissue. 


Two identical stationary electrodes placed about 2 cm 
from the subject would be unsatisfactory for hyperthermia 
therapy. There are two temperature peaks produced, one of 
which could over heat healthy tissue while a tumor was 
being treated with the other (Figure 4.7(a)). The solution 
would be to use one very large electrode (Figure 4.7(b)) or 
else asymmetrically spaced electrodes (Figure 4.7(c)). The 
latter case wasS examined in a manner similar to the 


Symmetric electrode case. 
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FIGURE 4.7 Hot spots produced in a dielectric by various electrode 


configurations. 
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The temperature distribution inside an agar cylinder 
irradiated with one electrode (2.0+0.1) cm, and the other 
(7,040.25) CM,. strOMe ses Surface (Figure 4.6) 15 shown in 
figure 4.9. There is only one large temperature peak 
produced with this configuration. For a shallow tumor this 
would be a Superior applicator to the symmetrically placed 
electrode type. There would be no danger of producing an 


unwanted hot spot on the far electrode side of the subject. 


The asymmetric electrodes produce a slightly lower 
peak temperature (0.2°C less) than the symmetric electrodes 
and a 0.9°C lower centre temperature. Far less energy is 
being deposited in the agar while little change is made in 
the peak temperature. Obviously the asymmetric electrodes 
are very efficient with regard to heating a shallow tumor 


and not imposing an unduly large heat load on the patient. 
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MEGAERA SIMULATIONS 


A considerable amount of labour was involved in 
producing each of figures 4.5 and 4.9. To expedite matters 
Megaera was used to simulate the temperature distributions 
in Tyei ice irradiated agar. Them effects Vote varying 
experimental parameters were quickly seen with computer 
Simulation. To validate the use of Megaera in place of 
actual experiments comparisons were made between measured 


and calculated temperature distributions. 


Results of symmetric electrode simulations made using 
Megaera are shown in figure 4.10. Figure 4.10 (a) and (c) 
are the maximum and minimum temperature limits calculated 
by Megaera when uncertainties in measured power, thermal 
and electrical conductivities, electrode spacing, 
temperature meaSurements, and volumetric heat capacity are 
taken into account. Figure 4.10 (b) is the actual measured 
temperature profile. The measured profile sits well within 
the maximum and minimum Limit si The temperature 
distribution is very sensitive to electrode spacing, power, 
and heat capacity. AS a reSult, the maximum and minimum 


temperature limits are widely spaced. 


A close approximation to the measured temperature 


profile is obtained from Megaera when the following 
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PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 


(a)-Maximum 


(c)-Minimum 


(b) -Measured 


28 


TEMPERATURE (°C) 


bo 


20 = 
DISTANCE (cm) 


FIGURE 4.10 Temperature profile measured in agar with high and low 
limits calculated by Megaera. 
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Parameters are used (Figures 4.11 and 4.12), 


ee 

© scary |- 20 9/7 

C4 e072 Cem: 

k=0247 eW/mcCc 

Power-12.6 W (see Appendix #2) 


Electrode distance from agar surface-1.85 cm 


These values are all within the measured bounds of 
uncertainty. Temperatures near the centre of the agar are 
calculated high by Megaera because the program assumes 
infinitely long electrodes. At the centre of the agar 
cylinder this assumption of no field spreading in the axial 
direction could amount to a 7% increase in temperature 
(Appendix #3). The error associated with using a two 
dimensional model decreases closer to the electrodes. This 
trend is evident in figure 4.10 and 4.11 and could account 
for the broad temperature peaks predicted by Megaera. In 
addition, the thermal con@uctivity is likely high resulting 


in low temperature gradients in the calculated data. 


Changing the position of the surface used by Megaera 
to calculate the current integral produces variations in 
Calculearedumbemperaturesmmot mUpatoLO. iC. lence a=. 0.. °C 
uncertainty was assumed for all temperatures calculated by 


Megaera. The+0.1°C uncertainty takes into account a 
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PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 


O -CALCULATED 
D -MEASURED 
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FIGURE 4.11 Temperature profile measured in agar and calculated by 


Megaera. 
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possible programming error not the error associated with 


using an approximate means of calculation. 


Within certain limits, as outlined above, Megaera is 
capable of accurately simulating heating in an agar 
Cylinder ~exposed) §to 1.f.) (32/75 MHz or lower) energy from 


Symmetric electrodes. 


Asymmetric electrode simulations were also made using 
Megaera. The results showing the maximum and minimum 
temperature limits are shown in figure 4.13 and 4.14 (a) 
andue cy. Figure 4.13 (b) shows the measured x-axis 
temperature profile and figure 4.14 (b) shows the measured 
y-axis profile. In figure 4.13 the measured values lie 
well within the simulated maximum and minimum limits. The 
Simulation does not show the very slight temperature 
increase in the distant electrode side of the agar. This 
could be the result of uSing a high thermal conductivity 
which would tend to smooth out temperature peaks. Also, 
the temperature peak was so small, practically within the 
bounds of experimental error, that FOoUunad-Oort “and 
approximation error could “have deleted it from the 


Simulation results. 


Along the y-axis (Figure 4.14) the measured 
temperature profile sits between the maximum and minimum 


calculated limits. The high measured readings at the 4.0 
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and 4.8 cm distances indicate that the thermal conductivity 


used in the simulation was too high. 


Figure» 4.15 (a) shows the measured asymmetric 
electrode temperature profile along the “x-axis. Figure 
4.15 (b) shows the analogous profile as _ predicted by 


Megaera with, 


© water bath = 2° S/m 
1.96 S/m 


° agar” 
Ge dnixi 06d / °C mn: 
k=0 547 W/m°G 
Power-13.94 W 


Near electrode distance-1.85 cm 


Far electrode distance-6.8 cm 


Dhesestwo  tenperatiire  rotites along the x-axis ane = very 
close to one another. Figure 4.16 also shows how close the 
measured and simulated temperature profiles are along the 
y-axis. The high centre temperatures in the simulation are 
likely due to the two dimensional approximation (Appendix 
#3) The steep temperature gradients in the measured plot 
and the lower peak temperature in the simulation indicate 
that the thermal conductivity measured for the agar is 
high. As with the symmetric electrode case, however, the 
plots are sufficiently close within experimental limits to 
show that Megaera can closely simulate this kind of heating 


(Figure 4.17). 
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SIMULATION OF MOVING ELECTRODE HEATING 


Applying the sinusoidal motion weighting function to 
the data obtained uSing static symmetric electrodes 
produced the temperature distribution in figure 4.18. The 
Simulated temperatures were compared to measured 
temperatures along the 0°, 45°, and 90° radial lines 
(Figure 4,19). Temperatures along the 0° line were too low 
(Figure 4.20) while the temperatures along the 45° (Figure 
ac Wwepand 9908 — (Figure 64.22) lines were too high. Water, 
resistance impeded electrode motion enough so that the 
speed variation represented a flattened sinusoid, or 


semi-circle, rather than a perfect sinusoid. 


A circular variation in electrode speed results in the 
teniperat Ure wdiStrabutrionein Figure «4.237 Along the 0° 
radial the temperatures are very close to those actually 
measured (Figure 4.20). The same holds true along the 45° 
(Figure Ase) and 90° “(Figure 4.22) radials. The 
assumpeton Otwa) 'Clirculam svarivativon Im electrode speedy ais 
a) good one “producing Sresultis very close lo the measured 


values. 


A flattened sinusoid can be approximated by averaging 
the circular and sinusoidal weighting functions. The 


temperature distribution which results is shown in Figure 


103 


ee ML nites 
lenges os tie 
eenit falour “Ot of ws 
Woh dha avev soth “Oaiatynote crude tegter: afify 
evuigit) ‘se vt oioge asigists (042 odd ofidy. (ot, 6 4 
sina pict vee Grae Wankt (7o.) “owpeh, S00 ben 
ats dale oe: dyvdn’ Métion <“osccale Geben! e2ea3 
se (Biseniie feaestets 6 Ge newe tga ood teisav - . 
Sone davliseg + Adds INO eerie lea - 
: _ 
$42.4) e3linn1 vee Whe ins at s@beatesy taluodts A 
“0 ef? yauls ie.¢ «wisi 9. eeltpsaineid stusaysgen a 
RSiayd 7s eeun! ~- yewds (39 sie pe esepiay sec tataer 7 : 


A ete wheis bux2 ablet) an : 96.9 erupt?) 6ehoneen 
oar visite: sS-8 ie) ‘he (15.8  etuBrgy 
BE. Beaue edd: ineale a4 FOrsDI: wwivetta” @ $4 tos sQRUaee 
i soGis bet esli.o> patedtese’, @6e Ooeg Lia 
eeu fev 


grigeis 2 te Ieniaoiggs 9G 420 ciseesle SGenesgals 
OH aha) fa4e8 gal she hit tohidbwels the -meteiels Aes 
Gigi? 1! ovotd ai adiveds idine eootumlagede 


104 


°01e 


OF ysnoiy} uopyjzerazea paeds Tepfosnuts Zuyos Aepun seporjodeTe dT13euMIAS 
SuTIFqIO YAM paqeay Aeapuf{Ao ze8e oprsuy uoFInqTAAsSTp sinjetedwsaz pazetnoyTe) gt’ qNITA 


= =» 


= 
> 


- ra 


i stecetar capes tore emmsiony oom Sees ee 


“~e 


105 


wo (S0°0*S6°T) 


eee | 


*(meTA do) AaputTTAd ae3e u~t soTdnodow1ey perequnu JO WUeUSDe Tg 61’ FUNDA 


wo (S0°0FS6'T) 


YUAGNITAD UVOV 


—— TVIdVa 906 


eare wrt) > Sage ce) (leer -ieesdt pwregee ts 2 eames 
hess S} 2e6efle- 
| Ve] “ 
4 aes 
6 


TEMPERATURE (°C) 


106 


PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 
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FIGURE 4.20 Temperature profile along on radial. 
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FIGURE 4.21 Temperature profile along 45 radial. 
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TEMPERATURE (°C) 


PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 


FIGURE 4.22 
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Temperature profile along 90 radial. 
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4.24, Along the 0° radial the temperatures are slightly 
low (Figure 4.20). Along the 45° (Figure 4.21) and 90° 


(Figure 4.22) radials the temperatures are slightly high. 


The preceding evidence “indicates that “the circular 
electrode speed variation most closely resembles the actual 
electrode speed. Within experimental limits the flattened 
Sinusoid weighting function also approximates the measured 


temperatures, though not as well as the circular function. 


Temperature patterns for moving asymmetric electrodes 
have also been calculated using the three averaging 
functions. These profiles are compared, as with the 
Symmetric electrodes, to temperatures measured along radial 
lines in the agar. These radial lines are at 0°, 45°, 90°, 
and 135° (Figure 4.25). Very little temperature variation 
eceurs over “ares centered on the 180° radial’ line so this 


data iS not included for comparison. 


Figure 4.26 shows the temperature distribution which 
would occur if the asymmetric electrodes travelled with a 
Sinusoidally varying speed. Along the 0° radial the 
predicted temperature profile is too low compared to the 
measured profile (Figure 4.27). As with the symmetric 
electrode case, this shows that the actual electrodes are 
travelling with a slower peak speed than the sinusoid 


assumption would indicate. Hence, the electrode speed 
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FIGURE 4.25 Placement of numbered thermocouples in agar cylinder 
(top view). 
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FIGURE 4.27 Temperature profiles along 0 radial. 
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variation is more closely approximated by a_ flattened 
Sinusoid. Along) tiemed5? (Figure ©4528) Vand s90° = Fagqure 
4.29) radials the temperatures are slightly higher than 
those measured. Along the 135° radial (Figure 4.30) the 
measured and calculated temperature profiles are 


essentially the same regardless of weighting factor used. 


The temperature distribution resulting from a circular 
variation in asymmetric electrode speed is shown in Figure 
4.31. This motion produces e temperature profile along the 
O*) radial” (Figure 4.27/)ewkich is slightly higher) than the 
measured profile. Along the 45°, 90°, and 135° radials 
(Figures 4.28-4.30) the temperature profiles are close to 


the measured profiles. 


The temperature Gistribution: resulting, frome sthe 
flattened sinusoid weighting function is shown in Figure 
Ot Siar Calemlated. protilessraiong the 70>, 452) “and aisoe 
radial lines-(Figures 4.27, 4.28, and-4.30)-are- very close 
to the measured ones for this type of motion. Along the 
90° radial (Figure 4.29) the calculated profile is high but 


within experimental error compared to the measured profile. 


The measured temperatures in Figures 4.27-4.30 are 
often found to be lower than the calculated temperatures. 
Round off errors in the computed and measured values could 


account for some of this. A comparison of the temperatures 
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PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 
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FIGURE 4.28 Temperature profiles along 45 radial. 
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FIGURE 4.29 Temperature profiles along 90 radial. 
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PLOT OF TEMPERATURE vs. DISTANCE IN AGAR 
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FIGURE 4.30 Temperature profiles along 135" radial, 
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measured at) «the centre ~of the agar phantom for the 
stationary and orbiting electrode cases revealed a 
difference in power absorption which could account for much 
of the discrepancy. During the moving asymmetric electrode 
measurements slightly less power was absorbed than during 
the stationary electrode meaSurements. This was the result 
of an inability to accurately maintain the power output to 
the phantom. The temperature difference was too small, 
probably about 0.05°C, to be picked up with the equipment 
used. Nontheless, a borderline measurement like 23.85°C 
Gizeam figure: 4:27) might have been | -23.,9°C swith. Slightly 
more power. These variations are well within the 


experimental error already outlined. 


Figures 4.27-4.30 show that the calculated temperature 
profiles are all equivalent to the measured profiles within 
experimental bounds of uncertainty. Consistent trends, 
however, show that the flattened sinusoid and circular 
weight functions produce better results than the sinusoidal 


funeELon. 


The averaging programs HTAVE.3, HTAVE.4, and HTAVE.5 
(See Appendix #1) use the sinusoidal, circular, and 
flattened Sinusoid welghting functions to calculate the 
moving electrode temperature distributions. The results, 
when compared to measured values, show that these programs 


can accurately simulate moving electrode temperatures by 
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using static electrode data. 


SIMULATED VARIATION OF APPLICATOR PARAMETERS 


The preceding sections establish the validity of using 
MEGAERA and the HTAVE programs to simulate stationary and 
moving electrode heating. This provides a means of 
observing the effect, on the temperature distribution, of 


altering various parameters. 


Figure 4.33 results from using twice the power as _ in 
Figure 4.12 with all other parameters the same. The change 


in centre temperature, 


where, 
T.-the temperature (°C) at the centre of the agar 
phantom. 
AT -the change in the centre temperature (°C). 
iSvabous 100% higher than win ebagure 4. 12). Similarly the 


change in peak temperature, 


where, 
T, the highest temperature (°C) cecurring am) the 
phantom. 
ATES change in peak temperature (°C) 


is about 100% higher. These results are expected in view 
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of the doubled input power. However, this linear relation 


is not obtained at higher power levels. 


Figure 4.34 results from using the same parameters as 
in Figure 4.12 except that the power is set four times 
higher to 50.4 watts. The temperature peaks occur at about 
the same positions in each case. AVE Be: 296% higher and 
Aalis Voees) \tonnighecm thane in ‘Figures. 451206 eihermal 
conduction into the water bath prevents Nagle and AT, from 


being linearly proportional to the power used. 


The effect of varying the water bath temperature is 
shown in Figure 4.35. The parameters used are the same as 
in Figure 4.34 except that the saline bath temperature is 
assumed to be 10°C. The temperature peaks are in the same 
PeSitions as with the 22. 1°C bath but AT, is 17.5% lower 
while AM is- only 1.8%. dower. This demonstrates the 
importance of having a cooling water bath. A substantial 
decrease in superficial temperatures can be achieved 


WEtnOULeStOniticantly altectring) deep temperatures. 


Figure 4.36 resultsetrom using the same parameters as 
in figure 4.306 ‘except “that. the “saline” bath) “electrical 
conductivity is set to be five times higher. The 
temperature changes Ale and NIE are both about 47% lower 
than in Figure 4.33. This demonstrates the importance of 


Usingeasbath oo bolus; swnien provides an electrical match 
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to the irradiated dielectric. Matched impedances result in 


maximum power transfer. 


Figure 4.37 shows again the effect of having a massive 
impedance mismatch. In this case the electrical 
conductivity of the bath is assumed to be about nineteen 
times lower than in Figure 4.33. The Ae is 75% lower and 


the Ae 80% lower than in the matched situation. 


Figure 4.38 demonstrates the effect of positioning the 
elecerodes! "35555 (cm, fromyche agar instead of 12685 cm. ~The 
AVE is reduced by 56% and the Ane by §41%. The smaller 
dO prin Al, implies that slightly improved deep heating may 


result from bnis Conf guranlvonmover that in figure «4.33. 


The change in temperature distribution produced by 
using progressively wider electrodes is illustrated in 
Figures 4.39 to 4.44. To produce these figures’ the 
parameters are slighty different from those used _ for 
Biaqures 4.2. 0r 4.335 SAS vamresult,. theyscannot be directly 
compared to those figures. Taken between themselves, 
however, Figures 4.39-4.44 serve to demonstrate the trend 


produced by widening electrodes. 


Figures 4.39 and 4.40 show the temperature pattern 
produced using 1 cm wide electrodes. Temperature peaks are 


higheeand close to the Strrace ms Ihescencrescemperature 1S 
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FIGURE 4.39 Calculated temperature distribution in medium using 
1 cm wide stationary symmetric electrodes. 
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FIGURE 4.40 Three dimensional plot of Figure 4.39. 
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low resulting in large thermal gradients. Figures 4.41 and 
4.42 show the pattern resulting from 2 cm wide electrodes. 
Temperature peaks occur slightly deeper than with the 1 cm 
electrodes. The ATpis 8.5% lower with the 2 cm electrodes 
and the AT. is 12% higher. A smoothing effect is seen to 
occur in the temperature distribution when wider electrodes 
are used. This effect is also demonstrated in Figures 4.43 
and 4.44. These simulations assume 4 cm wide electrodes. 
Again, the temperature peaks are broader and occur about 
0.5 cm deeper than in Figure 4.39. The Aa is 26% lower 
than with 1 cm wide electrodes and the AT. is 22% higher. 
For producing deep heating this levelling effect is most 


useful. 


Higher centre temperatures are produced by wide 
electrodes because the field between them is more uniform. 
This situation 1s roughly analogous to comparing the 
E-field between two point charges with that between two 
infanite Vine charges. The electric freld around a point 
charge varies with the inverse square of the distance from 
the charge. In the case of an infinite line charge the 
electric field is inversely proportional to the distance 
from it. Consequently the field between the lines does not 
vary as much as that between the two points. This is 
G@learly seen in Figures 4.3905 5to 4.44 as the electrodes 


become more like lines of charge rather than points. 
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FIGURE 4.41 Calculated temperature distribution in medium using 
2 cm wide stationary symmetric electrodes. 
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FIGURE 4.42 Three dimensional plot of Figure 4.41. 
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FIGURE 4.43 Calculated temperature distribution in medium using 
4 cm wide stationary symmetric electrodes. 
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FIGURE 4.44 Three dimensional plot of Figure 4.43. 
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A trend exists which indicates that, for stationary 
electrodes, deep heating is most effectively produced by 
uSing large electrodes. When using submerged electrodes 
Care must be taken to ensure that the r.f. power amplifier 
is not loaded down by the low impedance of large 
electrodes. A matching network can be used to protect the 
Eo oaeeanpliiier., However, there could be a problem with 
most of the power being dissipated in the matching network 
ana mot) the toads ‘This results in inefficiency but may not 


be important if sufficient power is available. 


Having established that the temperature peaks get 
smaller and the centre temperature higher as the electrodes 
get wider, it is logical to see if the centre temperature 


can be maximized in this way. 


Figure 4.45 was produced using the same parameters as 
Pugure 2333 except that the electrodes were set tos bes the 
Same width as the agar phantom. Changes in the temperature 
distribution were similar to those encountered in Figures 
4,39 to 4,44 with an important exception. Peak 
temperatures were Slightly higher (1.3%) than those 
encountered with Narrow electrodes. More power was 
dissipated between the wide electrodes, and less elsewhere 
in the saline. bath, than in the narrow electrode caSe. 
Further, the agar phantom essentially filled the space 


between the electrodes. This led to relatively high 
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temperatures throughout the agar. The temperature peaks 
were high because they occurred over 1 cm deeper in the 
agar than in Figure 4.33. (This mitigated’ the effect of the 
cooling bath on the peaks and allowed high temperatures to 
be obtained. HAY, was 13.4% higher than with the narrow 
electrodes. Stronger fields and smaller internal 
temperature gradients, as previously mentioned, were 


responsible for this result. 


To make the fields inside the agar even more uniform 
heating from 23.7 cm wide electrodes was simulated (Figure 
4.46). This produced one broad temperature peak in the 
centre of the phantom. The peak was quite low compared to 
the centre temperature in Figure 4.45 but higher than the 
centre temperature in Figure 4.33. The low temperature was 
the result of much of the power being dissipated in the 


Saline between the wide electrodes. 


For deep heating purposes the magnitude of the 
Pentrallyasproduced | “het spot’ a1s unimportant. Higher 
temperatures can be obtained simply by using more power. 
Whate1s important 1Ssthe fact that thedreatest temperature 
is produced deep within the agar and not near the surface. 
For stationary electrodes the oversized applicator produced 


the most desirable temperature pattern for deep heating. 
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ROTATING ELECTRODES 


To improve deep heating orbiting electrodes can be 
used. iiss setlLectw 1omecreanly Shown ink igqure «4.4 sere 
solid lines correspond to the temperature profiles measured 
along the x-axis (a) and the y-axis (b) for the stationary 
electrode case (Figures 4.5 and 4.6). The dashed lines (a) 
and (b) are the profiles measured along the x and y axes, 
respectively, for the orbiting electrode case (Figures 4.8 
and 4.9). In each case the electrodes are (1.9540.1) cm 
from the agar surface. The orbiting electrodes travel 
through an’ angle jot 240°, or an 80° arc, for these fiqures. 
This movement produces a AT» which is 30% lower than the 
Stationary electrode case with no change in the Nee 
Spreading the r.f. field exposure over the surface of the 
phantom produces lower, broader, slightly deeper peak 


temperatures than found with stationary electrodes. 


The effect produced by orbiting electrodes can be 
maximized by making the electrodes revolve completely and 
continuously around the phantom. Figure 4.47 is the 
Simulated result of using such revolving electrodes. Data 
from Figure 4.5 is used along with an averaging program to 
produce Figure 4.47. Temperature peaks are annular regions 
which have a ee 41% lower than the stationary electrode 


Case (Figure 4.5). ile ies the same aS with stationary 
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electrodes. Further, the annular peak is 2.2 cm deep 
inside the agar where as the stationary electrode peaks are 


only 0.9 cm deep. 


Figure 4.48, based on data from Figure 4.9, shows’ the 
expected temperature distribution from revolving asymmetric 
electrodes. The temperature peak is an annular region with 
a AT, 59% lower than in the stationary electrode case. 
This large drop occurs because the temperatures are 
averaged over one temperature peak instead of two as in 
Figure 4.47. There is no change in the centre temperature. 
The temperature peak occurs 1.5 cm deep within the agar 


instead of 1.1 cm with stationary electrodes (Figure 4.9). 


Applying the revolving electrode program to the data 
in Figure 4.45 results in Figure 4.49. The temperature 
peak is no longer an annular region but a broad peak in the 
centre of the agar phantom. Unlike Figure 4.46 this 
central temperature peak 1S quite substantial having 
relatively steep temperature gradients using the same 
amount of power. The revolving! 9.7 “cmp electrode 
configuration appears to be the best choice for producing 


deep heating in a 10 cm diameter agar cylinder. 


The model indicates that the best electrode 
Som roquracvon for producing deep hyperthermia uses 


electrodes which are the same width as the heated subject. 
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Also, the electrodes should be symmetrically placed 


relative to, and revolve continuously around, the subject. 


HYPERTHERMIA IN A LIVE ANIMAL 


A temperature profile obtained inside a live cat 
irradiated with orbiting symmetric electrodes is shown in 
Figure 4.50. Due to a lack of r.f. power the temperatures 
were quite low. The saline bath and matching network 
absorbed a significant amount of the energy from the 
electrodes. Blood flow may have reduced much of the 
temperature increase which would have’ occurred. Heating 
for longer than 20 minutes would probably have increased 


the core temperature but the cat could have awakened. 


The initial temperature in the cat was not actually 
uniform. The skin temperature was about 38.3°C due to the 


use of a heating pad during surgery and gradually fell to 


about 37.4°C deep inside the animal. Thisee Varcatton sin 
temperature likely influenced the final temperature 
GVstribucion, Time ~could = inot.| be takenehto allowr the 


animal's body temperature to reach a uniform steady state 
because of the danger of the anaesthetic wearing off. A 
cooler water bath was not used even though it would have 
controlled the subcutaneous fat temperature because the 
danger of producing hypothermia in the anaesthetized 
animal. The small temperature increase exhibited by the 
tied off kidney was probably due to surrounding blood 


vessels acting as heat sinks and a low initial 
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temperature. The temperature peak which occurred at the 
part of the kidney furthest from the electrode was found in 
two cats which were treated. A probable explanation was 


the lack of cooling bloodflow in this region. 


The peak temperature in) the) tat: (40,5°C)™ sror 
Stationary electrodes was estimated using the difference 
between moving and stationary electrode data in the agar 
(Figures 4.5 and 4.23). This was likely a conservative 


amount. 
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DISCUSSION 


EVALUATION OF ELECTRODE CONFIGURATIONS 


A figure of merit which could describe any 
hyperthermia applicator is the "target heating 


Speciiacity: ; 


AT 
> target 
n= 


T : 
non-target maximum 


where, 


-the temperature change produced at the 
target 
target site. 
AL ; -the greatest temperature 
non-target maximum 
change produced in the dielectric at a site other 
than the target site. 
Ideally the target heating specificity should be greater 
ehaneones INnis WOould™=imply that tne target, Or = tumor, 1S 
hotter than any other part of the irradiated body. Hence, 
therapeutically significant temperatures could be obtained 


at the tumor site without risk of damaging any other 


tissues. 


A value _ for f Of one. or slightly “less, “for an 
applicator Nearing  adamemcoulc eestill "Indlcarem ag, USerUud 


device. Vigorous blood flow in healthy tissue could 
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produce lower temperatures than implied by the agar or 
Simulations while sluggish circulation in a tumor could 
produce higher temperatures. AS a result, an applicator 
with an 74} in agar or Simulations could conceivably produce 


an>! ina live Subject. 


Plainly the symmetric electrode applicator which 
produced the temperature distribution in Figure 4.5 would 
give an ‘<I unless a tumor happened to exist at the site of 
each temperature peak. While the second temperature peak 
in healthy tissue would likely be smaller than that in the 
tumor, there would always by the risk of overheating 


healthy tissue. 


A more practical configuration would use asymmetric 
electrodes (Figure 4.9). The asymmetry could take the form 
of differently spaced electrodes, as used in this case, or 
different sized electrodes. Either case would produce one 
temperature peak aS in Figure 4.9. As long as the tumor 
waS positioned properly this type of applicator could be 
very useful. Deep heating would be beyond the ability of 
this type of applicator though. Some steering of the hot 
spot could be accomplished by varying the bath temperature, 
electrode size, electrode distance, and by moving the 
electrodes. However, the results would be difficult to 
predict, especially in an inhomogeneous human, and the hot 


spot would grow large and diffuse as depth increased. This 
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large hot spot would negate the need for steering. 


Very few applicators are capable of heating deep 
seated lung tumors. For these tumors an applicator with an 
ian 2) is requiredwewetaple 5.1 lists eraree values for 
most of the electrode configurations considered in this 
work. It should be noted that the (cenryre figures based on 


Megaera Simulations are likely to be slightly high due to 


the two dimensional nature of the program. 


The electrode configurations which produced an 
Lee > were the’ 23.7 temestatiuonary symmetricrelectrodes 
and the 9.7 cm revolving symmetric electrodes. In each 
case the temperature peak occurred in the centre of the 
agar phantom. Field Spreading could reduce these 
temperature peaks in a three dimension situation but the 
NeeNTRE value would probably change very little. The 
temperature peaks were broad and relatively flat. This 
would reduce the need for accurately aiming the hot spot. 
A tumor with poor blood flow situated anywhere in the 
temperature peak region would tend to heat up substantially 
more than its Surroundings. From this point of view a 
small, sharply peaked, hot spot would be undesirable 
because of ‘difficultyeinmaccurately positieoningsit inside 


an inhomogeneous human. 
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TABLE 5.1 


SUMMARY OF TARGET SPECIFICITIES (rn) FOR VARIOUS APPLICATORS 


n 


CENTRE 


0.49 


O83 


O55 


0.38 


ORTA! 


0.68 


Figure i 


4.5 


Description 


1 cm wide stationary symmetric electrodes, 
measured data 


1 cm wide stationary asymmetric electrodes, 
measured data 


1 cm wide stationary symmetric electrodes, 
simulated 


1 cm wide stationary asymmetric electrodes, 
simulated 


1 cm wide symmetric electrodes, sinusoidal 
speed variation orbiting +40 , simulated 


1 cm wide symmetric electrodes, circular speed 
variation orbiting +40 , simulated 


1 cm wide symmetric electrodes, flattened 
sinusoidal speed variation orbiting +40 , 
simulated 


l cm wide asymmetric electrodes, sinusoidal 
speed variation orbiting +40°, simulated 


1 cm wide asymmetric electrodes, circular 
speed variation orbiting +40 , simulated 


1 cm wide asymmetric electrodes, flattened 
sinusoidal speed variation orbiting +40 , 
simulated 


1 cm wide stationary symmetric electrodes, 
25.2 watts, simulated 


1 cm wide stationary symmetric electrodes, 
50.4 watts, simulated 


1 cm wide stationary symmetric electrodes, 
50.4 watts, simulated, 10 C bath 


1 cm wide stationary symmetric electrodes, 
i x 
25.2 watts, simulated, 5 Oy ath 
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TABLE 5.1 (continued) 


1 cm 
Zoee 


1 cm 
2S). 2 
from 


Description 


wide stationary symmetric electrodes, 
watts, simulated, o 19 
; fe bath / 


wide stationary symmetric electrodes, 
watts, simulated, electrodes 3.55 cm 
agar 


9.7 cm wide stationary symmetric electrodes, 


Ze 


em 
255 
Bane 


1 em 
25.2 
360 


watts, simulated 


cm wide stationary symmetric electrodes, 
watts, simulated 


wide rotating symmetric electrodes, 


watts, measured data, rotation through 


wide rotating asymmetric electrodes, 
watts, measured data, rotation through 


97 cm wide rotating symmetric electrodes, 


25.2 
360 


watts, simulated data, rotation through 
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Figure 4.49 shows a temperature distribution which is 
more desirable than that in Figure 4.46. The temperatures 
are much higher for the same amount of power which shows 
that the thinner electrodes are more efficient. More power 
1s dissipated in the agar with 9.7 cm electrodes than with 
23.7 cm electrodes. Using the revolving electrodes the 
temperature peak is not as flat as with the wider 
Stationary electrodes. This means a higher proportion of 
power 1s being deposited in the centre of the agar phantom 
than near the surface. The temperature peak could be made 
Steeper by using a cooler saline bath. This is borne out 
by companyison of POR ey: values for Figures 4.34 and 4.35. 
Use of a 10°C bath increased )..,,@¢ from 0.57 to 0.68. 
The need for having a bath or bolus for matching purposes 
aS “well Sdemonstrated anh Figures 4533 98 4536, "and' 4.37. 
Changing the electrical conductivity of the bath either 
substantially lowers YW ralayare (RUgurer Sem) mote nase) mele 
etiect ton 2t (Figure 4236) .eeln each of “Figures 4.369 and 
4.37 less power is deposited in the agar than when there is 
an electrical match between the bath and agar (Figure 
ARS ot Increasing the spacing between the electrodes also 
increases Ucenreg 2S in Figure 4.38 (i.e. Vee reae =0.75 
instead of 0.57 for Figure 4.34). The price paid here is 
less efficient use of power. Simply using electrodes the 
same width as the subject being heated (Figure 4.45) gives 
a much higher ere e Geer 09) and temperatures’ than 


changing electrode spacing. The improvement that results 
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from orbiting or revolving the electrodes around the 
phantom is clearly shown in Table 5.1. For the symmetric 
electrodes Y CENTRE is 0749 for Stationary electrodes 
(Figure 4.5) and 0.83 for the revolving electrodes (Figure 
4.47). Similarly for asymmetric electrodes (eee se 16) che: 
for the stationary case (Figure 4.9) and 0.8 for the 
revolving case (Figure 4.48). The increase in (ears 
because of orbiting electrodes is dependent on _ the 
variation in temperature around the agar phantom. In fact, 
because the temperature contours are very close to being 
Circular in Figure 4.46 moving electrodes would have 
virtually no effect on the temperature distribution. This 
result plus that with the widely spaced electrodes (Figure 
4.38) argueS against the utility of LeVeen's moving 
electrode hyperthermia applicator [38]. The results 
indicate that little benefit would be derived by moving the 


widely spaced electrodes about. 
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CLINICAL CONSIDERATIONS 


The advantages of uSing orbiting electrodes with 
surface cooling and dielectric matching are apparent for 
heating deep tumors. In addition, if parallel electrodes 
are used they should be as wide as the object being heated. 
Unlike agar the human body has a layer of poorly conducting 
fat just below the skin which tends to overheat in high 
power r.f. fields. Fortuitously orbiting electrodes reduce 
Surface, or near surface, temperatures the most. Combining 
orbiting electrodes with adequate Surtace™ "cooling, 
sufficiently wide electrodes, and enough power, high core 
temperatures could be obtained with tolerable fat 


LeMNpereacuLres. 


In a clinical setting more power could be made 
available to overcome the intrinsic inefficiency of this 
system. Large amounts of power can be absorbed by the 
coupling medium (saline) because it is lossy and acts as a 
heat sink. Glinveally thas power loss is not Signiticant 
as temperature distribution not efficiency 1s the primary 


consideration. 
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CONCLUSIONS 


The temperature distribution data clearly shows the 
advantages of using moving applicators to produce deep 
hyperthermia. Further, a ecoling bolus which is 
electrically matched to the subject and electrodes as wide 
as the subject are important. This should allow high 
temperatures to be obtained deeply without overheating 
Surface material. Clinically this means that deep tumors 
have a much better chance of being heated with the 
applicators described above than with the usual stationary 


applicavors.. 


Live animal experiments have not yet demonstrated deep 
internal heating. More r.f. power, or a smaller animal, 
coupled with a cooled saline bath should produce better 


temperature distributions. 


In summary, r.f. heating of agar has been simulated 
for both moving and stationary electrode cases. From this 
model an applicator using rotating electrodes the same 
width as the subject anda cooled, electrically matched 
bolus has been found to produce deep heating without high 


Superficial temperacunes. 
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APPENDIX 1. 


COMPUTER PROGRAMS TO SIMULATE R.F. HEATING WITH 


ROTATING ELECTRODES 


The following are printouts of the computer 
programs used to calculate orbiting electrode 
temperature Gistributions using stationary electrode 
data. These programs assume that the electrodes 
teavel back and forth through an bare son. plus-minus 
a0ee Throughout this arc of travel the electrodes 
vary in speed from a minimum at the ends of the arc to 
@ maximum at the centre. Ideally the variation in 
electrode speed would follow simple harmonic’ motion. 
This would mean a Sinusoidal variation in speed with 
time. For reasons discussed in the text this is not 
the case. A graph of electrode speed vs. time yields 
not a perfect sinusoid but a sinusoid with flattened 
peaks and troughs. In the extreme the variation in 
electrode speed was assumed to follow curve. which, 
over a range of 180°, was a half-circle instead of a 
Sine wave. An average of the ‘circular’ variation and 
the sinusoidal variation was used to give a ‘flattened 


Ginuscoldal. Vanlation. 


The program in Figure A.1.1 produced the 


temperature distribution for sinusoidally varying 
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electrode speed. 
Circular and 


respectively. 


Figures 


flattened 


Aioliize cand = AGtl eS 


Sinusoid speed 
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SINUSOIDAL ELECTRODE SPEED VARIATION 


REAL R(216),TH(216),T(216) ,TA(160) 
N=216 

READ CS 1OJeCR (i) aril i TC) 1=1,N) 
10 FORMAT(2F5.2,F6.2) 

M=5 

DO.20 121,160 

TA (I) =0704704*7(M) 405050074 (T(M-1) 47 (Ms 1))) 
*+0.06141*(T(M-2)+T(M+2) )+0.09409%(T(M-3) 
#+T(M+3))+0.2709%(T(M-4)+T(M+4) ) 
PRINT50,R(M),TH(M) ,TA(1) 

50 FORMAT(3F7.2) 

IF(M.EQ.23) GOTO 30 

IF(M.EQ.50) GOTO 30 

IF(M.EQ.77) GOTO 30 

IF(M.EQ.104) GOTO 30 

IP(M.BO. 131) GOTO 36 

IF(MIEO. 158)" GOTO) 30 

IF(M.EQ.185) GOTO 30 

IF(M.EQ.212) GOTO 40 

M=M?+ 1 

GO TO 20 

30 M=M+9 

20 CONTINUE 

40 CONTINUE 


STOP 
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END 


$ENDFILE 
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CIRCULAR ELECTRODE SPEED VARIATION 


REAL R(216)) THC 216), 7(216);TA( 160) 
N=216 

READ(S, 10) “(RUT ) Thi) T (1) t=17.N) 
LOO PORMAT(2ZP5e 2-6 Oz) 

M=5 

DOP 20 l= 1760 

TAC )=Cs 0S tts TOM) Om0eS2*.0T (M1). (Mat)) 
*+0.0905%(T(M-2)47(M+2)) +0. 1088*(T(M-3) 
*+T(M+3))+0.177%(T(M-4)+T(M+4) ) 
PRINT50,R(M),TH(M),TA(T) 

50 FORMAT (3F7.2) 

IF(M.EQ.23) GOTO 30 

IF(M.EQ.50) GOTO 30 

IF(M.EQ.77) GOTO 30 

IF(M.EQ.104) GOTO 30 

IF(M.EQ.131) GOTO 30 

IF(M.EQ.158) GOTO 30 

IF(M.EQ.185) GOTO 30 

IF(M.EQ.212) GOTO 40 

M=M+ 1 

GO TO 20 

30 M=M+9 

20 CONTINUE 

40 CONTINUE 
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$ENDFILE 
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FLATTENED-SINUSOIDAL ELECTRODE SPEED VARIATION 


REAL R216) TH (216), 02 lee TA 160) 
N=216 

READ (Sel 0 icra he (ts) el Cl)ad = iy, N) 
10 FORMAT. 2F 52,66. 2) 

M=5 

DO. 20) 1=49 160 
TA(T)=0706407*T(M) +0. 066635% (T(M—1) +71 (M=1)))) 
#70.0759555(T(M-2) 41 (M22)) +0. 101445407 (M—3) 
*+T(M+3))+0.22396%(T(M-4)+T(M+4) ) 
PRINT50,R(M),TH(M) , TA(T) 

50 FORMAT (GPa) 

LF(M.BO.23)) GOTO 2e 

IF(M.EQ.50) GOTO 30 

IF(M.EQ.77) GOTO 30 

IF(M.EQ.104) GOTO 30 

TECM2 BO, 131 5GOTO "30 

TP(M. £0,158)" GOTO 30 

IF(M.EQ.185) GOTO 30 

IF(M.EQ.212) GOTO 40 

M=M+ 1 

GOnTO.20 

30 M=M+9 

20 CONTINUE 

40 CONTINUE 


STOP 
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END 


$ENDFILE 


FIGURESA. I. 3 


APPENDIX 2. 
CALCULATION OF POWER ABSORBED BY LOAD 
Knowing the impedance across the matching network in 
series with the load and the impedance of the load alone 
the power deposited in the agar can be calculated. A 
Circuit representation for the r.f. heating apparatus used 
iMernisS Project 1S sshownwinerioGure A.25. 1 
From Figure A.2.1 
Z1+Z2=50/0° 0 
and by measurement 
dig = 09 al Oe O44 4 63.63) e: 
therefore, 


Z2,=(19.56-j83.63) 


The power dissipated across 2Z,+Z2 waS approximately 200 


Watts. Therefore 
I=P/V=200 W/100 VZ0°=2 Amps 70° 


where I is the current in the circuit, V is the voltage 
ACTOSS the two )l0adS = ©cqmeancumc>,. and 9b 92s) they power 
deposited in the two loads. So the power dissipated in 2, 


(the agar cylinder and saline solution) was 
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R.F. GENERATOR 


FIGURE Ao2.1” Equivalent circuit for hyperthermia applicator 


184 
P2=1?R2=121.76 Watts 


One tenth of the agar, saline bath, and electrodes were 
considered for all simulations so an applied power of 12.18 


Watts was used. 


APPENDIX 3. 


MAXIMUM TEMPERATURE ERROR RESULTING FROM 2-D SIMULATION 


Megaera was used to help determine the inaccuracy 
produced by applying a two dimensional program to a three 
dimensional world. To do this two heating simulations were 
made. In one case the electrodes and dielectric were 
closely confined by an insulating material (Figure A.3.1). 
The insulator kept current density uniform between the 
electrodes. In the second case the insulating boundary was 
far from the electrodes allowing the field to spread 
(Figure A.3.2). In each case power, electrode size, and 
dielectric properties were kept the same. The irradiated 
medium wasS given an extremely low thermal conductivity to 
minimize heat flow. Temperatures between the electrodes 
were found to vary by up to 7% between the two 
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INSULATING BOUNDARY 


DIELECTRIC | DIELECTRIC DIELECTRIC 


FIGURE A.3.1 Electrode - insulating boundary configuration for 
uniform current density 


INSULATING BOUNDARY 


— Ny 


DIELECTRIC 


FIGURE A.3.2 Electrode - insulating boundary configuration for 
non-uniform current density 


APPENDIX 4. 


PERMITTIVITY “AND CONDUCTIVITY” DERIVATIONS FOR DIELECTRIC 


MEASURING CELL[7] 


The equivalent circuit seen by the HP4815A vector 
impedance meter is shown in Figure A.4.1. Rp is the probe 
resistance, Lp is the probe inductance, Cp is the probe 
GCapecitance, Co 1S) the empty cell" “capacitance,  ©™ is the 
capacitance with the specimen inserted, and R is the 


measured resistance with the specimen inserted. Define 


C =C+cC#e#+#C 
t ee KO 


fousimplity the tollowing caleulativons. SRoeis typically 
0.542 so it can be ignored. Therefore the impedance of the 


Circuit 1s 


2 
2 R " au, R CS 
i T#R70"C, 7 Je me re, ¢ 
where is the angular frequency of the applied signal. 


From Figure A.4.2 the imaginary and real components of za 


can be specified by 


aR? C 
i = Ea ey ee 
zsind Os ItR2u Ch 1 
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FIGURE A.4.1 Equivalent circuit for dielectric measuring 
apparatus 


IM 


FIGURE A.4.2 Imaginary and real components of complex impedance 
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and 


where z=|Z| and is the phase angle of Z. From 2 


Tee Rowe ce ee R 
t zcosd 
O'G 
te ete 3 
Substituting Z.imto liqives 
zsind = On - WRC, zcose 4 


Substituting 3 into 4 gives 


L. 
zsind = > - (Rzcos¢ - z*cos*¢) ? 


So 


Gale ora 2a zsingd + z° 


zcoso 


The cell constant for the measuring chamber is defined as 


[7] 
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Substituting 6 into 5 gives 


Q(w*L * - 2wL zsind + 27) 
ele Sh cael Seen 


zcoso 


Or 


aan zcosd 
ca rica ae Wapicoe ay ay ee eR PRS EU ACT 
POCO, - 0 ee ez) 


Now from 1 and 2 


WRC 


OLS) = 2SinQe= eae = 
zsind GRO oF WRC, zcos> 


ONE 
WRC zcosd = On - zsind 


so 


WL -—- zsin 
5 > 


t  WRzcosd 


Substituting 5 into 8 gives 


wL - zsind 
(e} = 


t 2 


w(w*L - 2wL zsind + 2? 
D ) 
By definition 


Coa Crore C 
1 Pp fe) 
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so 


wL —- zsind 


var OUZE ag L, = 2uL zsin¢ ) 7 ee ui Sy) 2 


Substituting. 6 ante os gives 


WL - zsind 


a at =ats cm )}-— 


is wiz” + OL * = 2wL zsind) 8 
P p 


i= 


For frequencies below 5 MHz L<< 1 so 7 and 10 can be 


approximated by 
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